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ABSTRACT 


A  theory  of  gaseous  discharge  from  the  end  of  a  tube 
was  constructed  using  an  isentropic  model  with  account 
taker,  of  axial  gradients  in  the  state  variables.  On  the 
assumption  that  tne  flow  rates  from  such  a  tube  were  not 
appreciably  altered  by  the  presence  of  conventionally 
designed  muzzle  brakes,  formulas  for  the  forces  on  the 
brake  and  tube  were  obtained  for  brakes  of  various  design 

In  order  to  implement  the  computation  of  parameters 
associated  with  a  complete  brake  analysis,  a  digital 
computer  program  was  written  for  the  Royal  McBee,  LGP-30 
which  permits  one  to  perform  an  analysis  with  relative 
east.  This  program  is  ihcluded  in  the  report. 

A  comprehensive  bibliography  on  muzzle  brake  studies 
gun  induced  shock,  and  allied  fields  is  also  included  in 
this  report. 
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LISP  OK  SYMBOLS 


Due  to  the  large  number  of  parameters  considered,  the 
symbol!-  used  are  divided  into  three  classes:  those  ap¬ 
plicable  in  general  and  in  particular  to  the  period 
prior  to  gas  ejection,  those  concerned  primarily  with  the 
period  during  which  gas  is  discharged  from  the  muzzle 
(discharge  period),  and  lastly  those  concerned  with  the 
design  of  the  bra'te.  Any  symbol  not  listed  here  or  one 
having  two  meanings  is  defined  in  one  of  the  appendlcies, 
where  it  first  appears. 


PART  I 
GENERAL 

Symbol  Meaning  Dimension 


*  travel  of  projectile  from  beginning  inches 

of  rifling 

length  of  bore  from  beginning  of  inches 

rifling 

D  diameter  of  bore  inches 


rN 


A 

Acnct 

Vc 


radius  of  a  tube  of  cross  sectional  inches 

area  A;  also  the  radius  of  a  sonic 

nozzle 

O 

area  of  bore  plus  cross  sectional  in* 

area  of  grooves 

area  of  ammunition  band  in  contact  in^ 

with  tube  after  shear  deformation 

total  internal  volume  in^ 

volume  of  chamber  in^ 


Mp  mass  of  projectile 

i!c  mass  of  charge 

Mig  mass  of  igniter 

MT  -  Mc+Mig  total  propellant  mass 
Mr  mass  of  recoiling  parts 


lbm 

l^m 

lbm 

lbm 

lbm 


1 


Symbol 

Meaning 

Dimension 

T 

temperature  of  the  propellant  gases 

°R 

V 

velocity  (with  subscripts) 

ft /sec 

v . 

muzzle  velocity  of  the  projectile 

f t/sec 

V 

specific  volume  of  the  propellant 

ft3 

gases 

P  ^  1 

'  ~V~ 

density 

lbm 

ft3 

P  *  P(f  t) 

pressure  in  the  gas 

psia 

Pmax 

maximum  pressure 

psia 

t 

time 

secs 

t . 

time  at  which  base  of  projectile 
leaves  gun  tube 

secs 

B  ( t ) 

force  on  breech 

lb{ 

a 

speed  of  sound  In  gas 

ft/sec 

M  ~  X 
a 

mach  number  in  gas 

'I  -  cp/cv 

ratio  of  specific  heats  of  propellant 
gas 

universal  gas  constant  =»  1545 

lbf  ft 
mole  ' R 

m 

molecular  weight  of  gas 

lbm 

mole 

R  ■" 

gas  constant 

Jbf  ft 

m 

lb„  '  r 
m 

CP 

specific  heat  of  gas  at  constant 

Btu 

pressure 

Cv 

specific  heat  of  gas  at  constant 

Btu 

volume 

T^R 

J 

Joules  constant  -  777.5 

2 

lbt  ft 
Btu 

2 


Symbol 


Meani ng 


Dimension 


acceleration  due  to  gravity  ■  32.17 

heat  of  explosion  of  the  propellant 
charge 

heat  of  explosion  of  the  igniter 


ft/sec2 


V-  c 

^ie  if 

TTc 

L“).+ 


energy  involved  in  a  process 

kinetic  energy  of  the  projectile  at 
t  -  t0 

kinetic  energy  of  the  gas  at  t„ 

kinetic  energy  of  the  recoiling  parts 
at  to 

kinetic  energy  spent  in  engraving  the 
bore  and  in  forcing  the  projectile 
thru  at  t0 

energy  remaining  in  the  gas  at  t0 

total  heat  released  by  charge  and 
igniter  during  combustion 

heat  released  by  charge 
heat  released  by  igniter 

length  of  tube  of  uniform  cross 
section 

length  of  tube  of  uniform  cross 
section  having  an  internal  volume 
equal  to  Vr 

dimensionless  length 


inches 


inches 


ct  «  vt 


dimensionless  velocity 


/3  -  Vo 

aav  init 

t  "  R_  D  f  (n) 

Pb 


dimensionless  velocity 


dimensionless  pressure  at  t 


V  -  ^  (r\ )  dimensionless  specific  volume  at  t , 


T  -A,)  dimensionless  temperature  at  t0  ; 

also  used  as  an  angle  in  the  brake 
design 


degrees 
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PART  II 


Symbol 


r 

Pb  “  P 


B 
P 
M 

BT  ■  B 
H 

F„ 


cml 


aT  - 


r  -  aTt 

$  - 


LIST  OF  SYMBOLS 

PERTAINING  TO  THE  DISCHARGE  PERIOD 


Meaning 

axial  distance  measured  positively 
in  streamwise  direction 

radial  distance  from  axis 

pressure  at  the  breech  (written  in 
digital  computer  program  without  the 
subscript) 

pressure  at  the  muzzle 
breech  force 

cumulative  momentum  discharged 
cumulative  mass  discharged 
breech  impulse 

cumulative  stagnation  enthalpy  dis¬ 
charged  from  the  muzzle 

force  on  brake  in  axial  direction 
of  gun  tube 

force  on  brake  normal  to  axis  of  gun 
tube 

resultant  impulse  received  by  brake 
and  gun  tube  during  gas  discharge  in 
axial  direction 


Dimension 

inches 

inches 

psia 

psia 

lbf 

lbf  secs 
lbm 

lbf  secs 
Btu 

lbf 

lbf 

lbf  secs 


impulse  received  by  brake  and  gun  tube  lbf  secs 
in  the  normal  direction 


A  g  pb  intt  ;  dimensional  constant 
ab  init  Mr 

used  in  formulas  in  Appendix  III 
dimensionless  time 

£5 _  :  dimensionless  breech  pressure. 

Pb  init 

defined  in  Appendix  III 


secs 


-1 
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Symbol 


Meaning 


Dimension 


A  - 


Ar  - 


V  - 


»T  - 

1  - 

V 


e 


n 


A 


Br  I  dimensionless  cumulative 


A  pb  init 
breech  impulse. 


defined  in  Appendix 


III 


dimensionless  momentum 


A  ^b  init 

rate  of  discharge,  defined  in 
Appendix  III 


P  g 


dimensionless  cumulative 


*T*binit 

momentum  discharged  from  muzzle, 
defined  in  Appendix  III 

-i -I 


Mt  06 


8  A  Pb  init 


^dimensionless 


init 

mass  rate*  of  discharge,  defined  in 
Appendix  III 

M  ;  dimensionless  cumulative  mass 

discharged,  defined  in  Appendix  III 
• 

H _  ;  dimensionless  stagnation 

£g 

enthalpy  rate  of  discharge,  defined 
in  Appendix  III 


H _  ;  dimensionless  cumulative 

*r 

stagnation  enthalpy  discharged, 
defined  in  Appendix  III 

expansion  ratio,  i.e.,  ratio  of  cross 
sectional  area  of  gas  flow  at  any 
axial  station  to  area  of  muzzle 

speed-up  factor,  i.e.,  ratio  of  one¬ 
dimensional  velocity  of  gas  at  any 
station  to  velocity  of  gas  at  the 
muzzle 

empirical  divergence  angle  of  the 
flow  with  respect  to  the  axis  at 
entrance  to  a  control  volume  surround¬ 
ing  the  baffle  (Appendix  V) 


dimensionless 


radians 


5 


to  I 


vinbol 


Meaning 


Dimension 


N 

r__  1 
rN  rH 


dimensionless  axlai  or  stream- 
wise  distance 

dimensionless  radial  or  lateral 
distance 

coefficients  used  in  a  cubic  fit 
of  the  shock  boundary  of  the  jet 
for  cTm  <2;  Vm  = 

2  3 

"h  m  "  a+bl«fl'®  +b SCTm 


oij  .  /3 


<*y  /6’ij 


coefficients  used  in  a  quadratic 
fit  of  the  shock  boundary  of  the 
jet  for  CTm  >2; 

7\  m  '  a'*bl  <Tm  +b2  CT m 

coefficients  used  in  a  cubic  fit 
of  shock  boundary  coefficients: 
a  ~  a(p./p5o  ), 
etc, 

coefficients  used  in  a  cubic  fit 
of  shock  boundary  coefficients: 
a  -  a  *  (P0/P«p  ).»  b{  ”  (Pc/Po<>)’ 

and  bg  -  b^Cp^/p^) 
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PART  III 


LIST  OF  SYMBOLS 
PERTAINING  TO  BRAKE  DESIGN 


Symbol 


Meaning 


Dimension 


z 


y 


ol 


J 


i 


axial  direction  ot  brake  measured 
positively  away  irom  breech 

direction  normal  to  axis  ol  brake 
measured  positively  upwards  when 
gun  tube  is  horizontal 

angle  thru  which  gas  is  deflected 
with  respect  to  the  z-  direction 

angle  thru  which  gas  is  deflected 
with  respect  to  the  y-  direction 
in  a  plane  norma]  to  the  z-  axis. 
See  figure 4.2  in  appendix  IV. 


f  i-0, 1, 2  closed-baffle  port  areas,  used 
J  ’  with  subscripts  to  indicate  the 

baffle  referred  to  in  a  multi¬ 
baffle  brake.  See  figure  4.1 
in  appendix  IV. 

rj  ratio  of  the  mass  flow  rate  cap¬ 

tured  by  the  j  th  baffle  to  the 
mass  flow  rate  oT  gas  from  the 
muzzle 


CO 


ratio  of  the  mass  flow  rate  of 
gas  passing  the  .1  th  baffle  un¬ 
deflected  to  the  mass  flow  rate 
from  the  muzzle 

ratio  of  axial  brake  force  to 
momentum  rate  of  discharge,  P 


ratio  of  resultant  momentum  rate 
of  discharge  in  the  z-direction 
thru  a  control  volume  enclosing 
only  the  brake  orifices  to  the 
momentum  rate  of  discharge  from 
the  muzzle 

ratio  of  normal  brake  force  to 
momentum  rate  of  discharge 


inches 


inches 


radians 


radians 
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■Symbol 


Meaning 


Dimension 


v  ~  F7  momentum  index,,  i.e.,  the  ratio 

g —  of  the  axial  brake  force  to  the 

breech  force 


empirical  correction  factors 


1,2S3  angles  associated  with  free- 
periphery  brake  design.  See 
appendix  VI. 


radians 


zm  distance  from  muzzle  to  pro-  inches 

jection  of  rim  of  open  baffle 
on  z-  axis.  See  figure  5.1b 


in  appendix  V. 

ym 

radial  distance  at  zm  from  axis 
to  shock  envelope 

inches 

i 

inner  radius  in 
See  figure  5.1b 

open  brake  design, 
in  appendix  V. 

inches 

Q 

outer  radius  in 
See  figure  5.1b 

open  brake  design, 
in  appendix  V. 

inches 
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SUBSCRIPTS 


refers  to  conditions  at  projectile  ejection 

refers  to  ambient  conditions 

pertains  to  the  gas 

pertains  to  the  projectile 

pertains  to  the  recoiling  parts 

pertains  to  the  igniter 

pertains  to  the  charge 

refers  to  value  of  variable  at  the  breech 
refers  to  an  effective  value 
refers  to  an  average  value 
refers  to  a  maximum  value 

refers  to  initial  conditions  for  gas  ejection 

refers  to  a  cumulative  or  integrated  value 

dummy  index;  takes  integral  values  from  0  to  m 

dummy  index  representing  the  baffle  number; 
take  integral  values  from  l  to  N 

indicates  coefficient  in  power  series,  baffle 
number,  and,  in  general,  distinguishes 
membei s  <>!  a  set 


SUPERSCRIPTS 


refers  to  a  stagnation  value 

indicates  first  deviative  of  variable 
with  respect  to  time 

indicates  value  of  variable  at  time  gases 
become  sonic  at  the  muzzle 

indicates  a  mean  value 
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INTRODUCTION 


The  primary  purpose  of  this  study  was  to  establish  a 
reliable  and  convenient  means  for  analysing  the  internal 
ballistics  of  an  artillery  piece  during  the  gas  discharge 
period  and  for  evaluating  the  performance  of  an  associated 
muzzle  brake  of  one  of  a  number  of  different  configura¬ 
tions.  Several  standard  brake  designs  and  some  novel 
designs  were  considered  and  are  examined  separately  later 
in  this  report. 

A  secondary  motive  for  this  study  was  to  establish 
relations  for  computing  parameters  useful  in  other  studies 
such  as  shock  establishment  and  decay. 

The  method  of  approach  to  the  determination  of  in¬ 
ternal  ballistic  parameters  during  the  gas  discharge 
period  involves  several  assumptions: 

X)  that,  at  the  advent  of  discharge,  all  of  the 
propellant  is  burnt; 

2)  that  the  heat  released  by  the  reaction  may  be 
computed  by  taking  the  sum  of  the  product  of  the 
heat  of  explosion  of  the  igniter  and  its  mass 

t  d  that  of  the  propellant  and  its  mass; 

3) '  that  the  amount  of  gas  discharged  after  projectile 

ejection  necessary  to  make  the  gas  velocity  at  the 
muzzle  sonic  is  negligible; 

4)  that  the  gas  discharge  proceeds  nearly  isen- 
tropically ; 

5)  and  that  the  presence  of  the  muzzle  brake  does 
not  appreciably  affect  the  discharge  rate  from 
the  muzzle. 

To  obtain  accurate  initial  conditions  for  the  gas  dis¬ 
charge  period  the  following  method  was  used.  An  analysis 
of  energy  disposition  prior  to  gas  discharge  was  made 
(Appendix  I).  As  a  result  of  this  analysis  one  can  compul- 
the  total  thermal  energy  remaining  in  the  gas  at  the  advent 
of  gas  discharge.  From  the  latter  value,  one  can  compute 
the  average  temperature  and  pressure  at  the  start  of  dis¬ 
charge.  By  making  certain  assumptions  regarding  momentum 
transport  within  the  tube  prior  to  gas  discharge,  one  can 
predict  the  relationship  that  exists  between  the  average 
values  of  the  temperature,  pressure,  and  specific  volume 
and  the  values  of  these  parameters  at  the  breech.  This 
derivation  is  made  in  appendix  II.  These  values,  then, 
are  used  as  initial  conditions  for  the  nearly  isentropic 
flow  from  the  muzzle 


11 


With  the  initial  breech  pressure  and  specific  volume 
given,  one  assumes  that  the  gas  rapidly  becomes  sonic  at 
the  muzzle.  By  applying  isentropic  flow  theory,  one  can 
obtain  expressions  for  the  breech  pressure  and  breech 
force,  the  momentum  rate  of  discharge,  the  mass  rate  of 
discharge,  the  stagnation  enthalpy  rate  of  discharge,  and 
integrals  of  all  of  these.  See  Appendix  III.  It  was 
found  that  convenient  non-dimensional  expressions  for 
these  parameters  could  be  obtained  in  terms  of  a  dimension¬ 
less  time  7"  .  Evaluation  of  these  expressions  provided  a 
set  of  tables  suitable  for  use  with  dimensional  constants 
in  obtaining  the  dimensional  parameters  relative  to  the 
discharge  period.  A  computer  program  which  takes  certain 
tliermodynamic  data  and  computes  the  dimensional  parameters 
was  given  the  acronym  "DIPARDIP." 

DIPARDIP  also  requires  muzzle  brake  data  and  furnishes 
loads  and  impulses  on  the  brake  as  well  as  the  effective 
momentum  index.  Use  of  this  program  and  some  typical  re¬ 
sults  are  described  in  later  sections  of  the  report.  A 
special  feature  of  this  program  is  the  large  number  of 
different  brake  designs  which  can  be  analysed. 

In  order  to  clarify  what  is  meant  by  the  terms  "closed 
brake"  and  "open  brake"  designs,  the  reader  is  referred  to 
figures  0.1  thru  0.4  .  Appendix  IV  treats  the  analysis  of 
the  closed  brake  design.  Appendix  V  treats  open  brake 
design,  and  the  special  equations  for  the  free  periphery 
design  are  derived  in  Appendix  VI. 

In  the  next  section  all  of  the  equations  used  by 
DIPARDIP  are  presented  together  with  tables,  graphs  of 
tables,  and  graphs  of  some  of  the  equations.  These  graphs 
and  tables  could  be  used  for  a  hand  calculation  of  the 
equations  in  lieu  of  DIPARDIP.  For  the  derivation  of  any 
of  the  equations,  the  reader  is  referred  to  the  appendicies. 
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-070-2150-U?  I  ROCK  ISLAND  ARSENAL  ORDNANCE  CORPS  I  3  January  1962 
HOWITZER,  LIGHT,  TOWED,  105mm,  XM102.  UBU  MUZZLE  BRAKE. 
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Figure  0,2 

A  TRIPLE-BAFFLE, 
CLOSED  SYMMETRIC 
BRAKE 
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EQUATIONS 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 


0.7 

0.8 

0.9 

0.10 

0.11 

0.12 

0.13 

0  li 
0.15 


^ig  =  eig  Mig 
ft c  "  ec  «c 

-ft  -  ftH 

KEP  ■  M_ 

P  0 

2  g  J 

KEG  «  Mt  Vo  (const) 

2  gT 

const.  - 


+ 

2 

'  Vr  1 

I  l+vr  1  ,  where 

V  v 

(vj 

Vo  1 

i  M 

vr  -  Mt 

v7  sr 

+  Mp 

mt 

+  Mr 

2 


XER  -  Mr  v* 

2  g  -J 

KEB  -  .18  Acnct  (KEP) 

- s - 

Eg  - 14  -  2(kinetic  energies) 

Tav  init  *  Tg  ’  gi iJllH  1 

R  Mt 

Prv  init  *  Per  “  12  E  C'j  -1)  J 

Vr  1/2 

aav  init  (Tf  g  E  Tav 
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aav  init 


^  av 

^  av 


1.0495200  - 
.97960297  + 


.25021815^3  -  . 061082024 
.10274869^6  +  .19109947^3  2 
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0.16 

^av  “  <f 

0.17 

Pb  init 

0 .  ^  8 

Vb  init  : 

0  jl9 

ab  init 

0.20 

C[T  **  A 

*b 

0.21 

t  -  T 

oCr 

0.22 

Pb  =  P  “ 

0.23 

B  -  A  pb 

0.24 

a)  Bt  »  A  p, 

b)  Bt  -  ab 

0.25 

P  -  A  pb 

0.26 

p  “  ab  i: 

av 


'If. 


av/  init 
-) 

^  av/  init 

aav  init 
<l/.v>  1/2 


3b  init 

r,T'¥_ 


<XT  T 

Lnit  mt  ^ 

g  T 


0.27 


g 

M  -  I g  A  ^ b init  \  y 
ab  <nit  j 


0.28  M  -  M  T  VT 

0.29  H  -  Eg  C*-r  r\ 

0.30  H  -  Eg  t^T 
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EQUATIONS  FOR  SYMMETRICAL 

AND  ASYMMETRICAL  CLOSED  BRAKES 

0.31  a) 

-  J S0j  +  Sij  ’ 
\  S0j-1  / 

where  b) 

6I0 

£  1 

c) 

s00 

s  A  and  j  is 

the  number  oi  the  baffle 

0.32 

£ 

2j 

Slj  J 

0.33 

nlj 

»  n(€1;J)n 

Function  is  tabulated. 

0.34 

n2j 

-  n(£2j)J 

0.35 

Vlj 

'  ao  nlj 

0.36 

V2j 

“  ao  n2 j 

0.37 

«2j 

*  rj  Mo 

0.38  a) 

r  , 

./  Slj  ) 

Kj-l  ,  where 

j 

lS0j  +  5lj  / 

b) 

Ko 

-  1 

0.39 

KJ 

-/  S0J  ' 

)  Kj-1 

\50j  ’  Slj  / 

/  J 

0.40 

Xr 

-  £  Q12 j  rj  cos  0<  jj  +  (1-2 

f2 

^r 

-  ^aps  -  n1N  (1-bps  jj/f2 

0.41  a) 

aps 

“  Jj  r j  n2j  cos  j 

b) 

bps 

N 

j“l  J 

c) 

cps 

-ji  r  j  n2 J  cos  /3  j 

d) 

dps 

-  I  r,  rilf 
J-l  J  lj 

j )  or 
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0.42 

"  -jll  rJ  "2J 

c°s  p  or 

fl 

C*J  •=  cps 

m- 

0.43 

^  -8  rj  (nl. 

i  ~  »2 j  cos  1 ) 

J**  J* 

t'o 

Cx  /»  -  upb  -  aps 
fo 

0.44 

>-  i-> 

0.45 

r2  -  ># 

0.46 

Fz  -  |?  B 

0.47 

Fy  -  COP 

0.48 

°*  ®T  “  ^  P 
Iz  -  (1  -  |?)bt 

or 

0.49 

ly  »  OJP 

EQUATIONS  FOR  SINGLE -BAFFLE 
SYMMETRICAL  OPEN  BRAKES 

0.50 

b  -  n2i  rx 

cosA(l  -  cos  c*  1) 

where  cosA  =  1 

0.51 

rl  -  1  -  I2  , 

f°r  ym  <  Q 

or  ‘I  ”  (1  -  +  (jJ-  q2  _  j2j  ^ 

ir  n 

y2 

'm 

0.52  1121  “  n 21^ 

0,53  ^21  “  1  +  3^i  Xm 


for  ym  >  Q 
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0.59 


0.54a 

0.54b 

=  ym/rH 
r..  =  (A/tt)1/2 

0 . 55 

<T  _  =  Zm 
m  'N 

0.55  a) 

7s_  =  a  +  b-i  CT  •+ 
''m  J-  m 

b) 

"X  =  a*  +  b’  <T 
Am  1  m 

0.57 

po  "  pb 

p~  TrTyyp^- 

0.58 

a  =  +  ft10  / 

2  +  *>3  <T  l  ,  for  ^ 

-  m 


m 


m 


bl  “  °*1  +  /^  11  /Ea.)  +^2l  /P«-\2  +  /^3i  /Pi) 

'  W  ;  \pJ  '  IrJ 

b2  -  ^2  +  /^12  f  Po)  +  & 22  /poj  2  +  /^32  /Po\ 

Viw  va »/  "  Ir»/ 

*/*“(£)  +  Z^33  (fe)  2  *  .^33  f^) 

*^30(t) 


3 

a' 


b’  -  o<;  + 


M 


;  Vii(fe)  V»(g)2  (g) 

i  +  ^2(ir)  +  <^2(fej2  V“(£) 


0.60  Iz  -Z  (1-  l^i)  (BTi-BT  i.j)  or 

lz  “  aTN“f.1^i(B'rl”BTi-l)  ’ 
where  H^q  =  0 

0,61  ^  eff  l?i(BTi-BTi_i> 

N 


0.62 


Iz  *  &rN  (l-)?eff) 
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0.63  00  =>  0  ,  for  a  symmetrical  brake 

0.64  7\r  =  n21  r1  cos  0<  j  +  n21  (l-r1) 

EQUATION  FOR  FREK-PERIPHERY  BRAKE 

0.65  60  -  1^  S  n2j  r j [cos i/i j-cos 2/3 j +sin  (ff/2— 

7T  j  L  J 
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.91138 

.20 

.83134 

.30 

.75899 

.40 

.69352 

.50 

.63422 

.70 

.53169 

1.00 

.41052 

1.30 

.  31910 

1.60 

.24962 

2.00 

.18162 

2.50 

.12377 

3.10 

.079635 

3.70 

.052227 

4.50 

. 030576 

6.10 

.011365 

7.50 

.0051568 

M- 

V 

.62049 

.68128 

. 56600 

.62736 

.51674 

.57815 

.47216 

.53321 

.43179 

.49213 

.36199 

.42013 

.27949 

.33314 

.21725 

.26575 

. 16995 

.21321 

. 12365 

. 16029 

.084268 

. 11364 

.054217 

.076510 

.035557 

.052405 

.020817 

.032418 

. 0077378 

.013343 

.0035109 

.0065672 

Table  0.  1 
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.095494 

.065014 

.071039 

.20 

. 18256 

.  12429 

.  13643 
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. 26202 

.  17839 

. 19667 

.40 

.33459 

.22780 

.25220 
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.40093 

.27297 

. 30344 

.70 

.51719 

.35211 

.39445 
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.44775 

. 50688 

1.30 

.76646 

. 52183 

. 59629 

1.60 
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.57958 

.66781 
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2.50 
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3.10 
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1.11096 

.75637 

.90408 
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6.10 
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.97143 

7.50 
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.76055 

. 16866 
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.68788 
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.62273 
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. 56427 

.36578 

.70 

.46454 

.46831 

i .  00 

. 34025 

. 58948 

1.30 

.26452 

.68090 

1.60 

.20176 

.75038 

2.00 

. 14207 

.81838 

2.50 

.093072 

.87622 

3.10 

.057223 

.92037 

3.75 

.035934 

.94777 

4.50 

.019909 

.96942 

6.10 

. 0066835 

.98863 

7.50 
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. 99484 
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. 92324 

5047 

.90800 

5561 

. 89197 

6077 

.87532 

6595 

. 85821 

7116 

.84079 

7638 

.82321 

8163 

.80557 

^  av 

av 

1.06639 

.98453 

1.08088 

.98144 

1.09665 

.97818 

1.11364 

.97479 

1.13177 

.97129 

1.15097 

.96772 

1.17117 

.96411 

1.19230 

.98048 
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1.5000 
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1.6000 
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1.8000 
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1 . 9000 
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2.5000 
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Table 


USE  OF  DIPARDIP 


This  program  was  designed  to  be  used  with  an  inter¬ 
pretive  routine  for  the  LGP-30  computer  called  DICTATOR. 
Consequently  the  conventions  established  for  that  routine 
will  be  followed  here. 

As  explained  in  the  INTRODUCTION,  DIPARDIP  will 
compute  the  various  parameters  of  interest  in  connection 
with  a  particular  weapon-brake  system  during  the  gas  dis¬ 
charge  period.  These  parameters  are  labeled  when  printed. 
A  sample  format  with  explanations  is  shown  on  the  next 
page. 
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Bop 

p  p  p 

09 

4) 

O  41  3 

49  3  it  3 

1  0 

ft  a  ^ 

41  ft 

P  09  P  ft 

1  p 

■M  -H  O 

m  u  a 

t>  X  B 

1  u 

•P  P 

43  P 

*  Vi  at  P 

w 

s/ 

w 

N  « 

1  Vi 

p 

.C< 

rH 

w 

n 
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To  obtain  these  results,  it  is  necessary  to  prepare 
a  different  type  or  data  tape  fcr  each  type  of  brake  con¬ 
sidered,  They  are  classed  as  follows: 

Type  1.  Closed.,  asymnetric  brakes 

2  Closed,  symmetric  braxes 

3  Open.  symmetric  brakes 

4,  Clcsedt  free -periphery  brakes 

For  each  of  the  types>  the  data  pertaining  to  the 
woapcn  and  propellant  charge  must  be  entered  on  the  data 
tape  as  follows: 


20000100' 

R 

ec 

e~g 


Acnci 

Yr 

Mq 

M 


Mr 


ig 


exit  ‘ 


The  data  tapes  for  each  type  may  be  completed  as 
shown  below,  (Refer  to  the  appropriate  appendix  for  a 
drawing  of  the  brake  ) 


Type  1 


20000070’ 


S0L  sii 

S21 

P  1  (by  rows) 

S02 

Sl2 

S22 

G<  2 

/S2 

S03 

S13 

S23 

C<3 

h 

exit'  ** 

*  All  numbers  must  be  written  in  the  "v50"  system  which 
is  used  with  the  DICTATOR  interpreter, 

**  The  "exit'  "  must  be  typed  after  j3  ^  ^or  J“1  or  2  for  a 
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Type  2 


The  same  data  tape  format  is  followed  for  this  type 
as  for  type  1,  with  the  exception  that  &  i  is  assigned 
the  value  1.5707,  jr/2  .  '  J 

Type  3 


20000090 • 

4  * 

I 

Q 

zm 

*1 

exit' 


Type  4 


In  addition  to  the  data  prepared  for  a  type  1  brake, 
the  following  data  should  be  entered: 


20000052' 


v  11 

v 21 

^31 

^12 

^22 

32 

*13 

■*43 

•^3 

exit’  **  . 

(by  rows) 


The  value  of  j  used  for  a  type  1  brake  is  not  used 
for  a  type  4.  Substitution  of  0  for  this  value  is, 
therefore,  recommended  in  preparing  a  data  tape  for  a 
type  4  brake. 


one  or  two  baffle  brake,  respectively,  i.e.,  after  the 
last  value. 

For  a  two-  or  three-baffle  brake  the  following 
must  also  be  typed: 

10001512' 

8000000(j-l) ’  . 

*  op  clt 

**  The  "exit'  "  should  be  placed  after  the  last  . 
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For  each  of  the  four  types  of  brake  designs  con¬ 
sidered,  there  is  an  initial  program  location.  Corre¬ 
sponding  to  each  initial  location,  there  is  a  start 
compute  instruction.  These  are: 


for  type  1. 

2. 

3. 

4, 


3000150^' 

30001500' 

30001000' 

30001505’ 


The  program  uses  all  breakpoints  except  8.  Therefore, 
these  must  be  depressed  for  uninterrupted  operation.  If, 
and  only  if,  the  high  speed  punch  is  to  be  UBed,  DICTATOR 
requires  that  break  point  8  be  depressed  also. 

When  the  break  points  are  sot  depressed,  the  com¬ 
puter  will  stop  —  at  conditional  stop  4  during  the 
computation  of  the  momentum  index  for  closed  brakes  and 
prior  to  printing  the  parameters  for  the  present  value 
of  time,  t  —  at  conditional  stop  16  prior  to  computing 
the  parameters  for  the  next  instant  in  time  —  and  at 
stop  32  prior  to  initializing  non-dimensional  parameter 
locations  for  type  1  brakes. 

When  used  with  the  high  speed  punch,  the  computer 
takes  about  twenty  minutes  to  complete  the  computations 
for  one  weapon-brake  system.. 
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Dl  PAR  Dl  P 


\ar*' 


iCrK.’o 


50 


DETAIL  LOGIC  DIAGRAM 
FOR  OPEN  BRAKE 
COMPUTATIONS 
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STORAGE  ALLOCATIONS 


Locations 

0012-0036 

0038-0050 

0052-0068 

0070-0098 

0090-0098 

0100-0128 

0130-0142 

0144-0174 

0176-0206 

0208-0238 

0240-027C 

0272-0302 

0304-0334 

0336-0366 

0368-0398 

0400-0498 

0500-0558 

0560 

0562-0648 

0650-0680 

0682-0738 

0740-0796 


General 


Contents 


constants 


open 

J. 


11 


^33  ,  for  free-periphery  brakes 


Sqi  . • • /^3  >  for  closed  brakes 
I,  Q,  s 


cK  ,  for  open  brakes 


interior  ballistic  data 
constants 
T"  table 
table 
jjl  table 
y  table 
<§T  table 
^Lr  table 
yr  table 


7r 


table 


working  storage 
print  out 
P/B  «  A/ $ 

-  .6808226  for  T  -  1.26 


open 

table 
n  table 
€  table 
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Locations 

Contents 

0798-0820 

<*6  •  •  •  /?32 

Set  I 

0822-0844 

<*0  '  ‘  '/^32 

Set  II 

0846-0896 

open 

0898-0928 

*0  "  *  ^33 

Set  I 

0930-0960 

<*0  ‘  '  ‘  /^33 

Set  II 

0962-0998 

open 

1000-1632 

program 

Print  Out 

Location 

Content 

0500 

t,  secs 

0502 

pb,  psia 

0504 

B,  lbf 

0506 

P,  lbt 

0508 

M,  lbm/sec 

0510 

Bj,  lbf  secs 

0512 

P,  lbf  secs 

0514 

M,  lb 

a 

0516 

H,  Btu 

0518 

*7/ ,  Btu  and 

H,  Btu/sec 

0520 

KEP,  Btu 

0522 

KEG,  Btu 

0524 

KER,  Btu 

0526 

KEB,  Btu 
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Location 

Content 

0528 

Eg,  Btu 

0530 

Tav  init. 

0532 

pav  init. 

0534 

1-  b 

0536 

b 

0538 

0540 

CO 

0542 

0544 

f 

0546 

Xz 

0548 

Fz 

0550 

xy 

0552 

Fy 

0554 

open 

0556 

0558 

,  b  eff 
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Constants  and  Internal 
Ballistic  Data 


Location 

Content 

0012 

0 

0014 

14.7,  Poo 

0016 

3.141593,  v 

0018 

40 

0020 

1.570796,  jr, 

0022 

1.0495200 

0024 

-.25021815 

0026 

-.061032024 

0028 

.97960297 

0030 

. 10274869 

0032 

. 19109947 

0034 

156 

C036 

1.33,  f3 

0100 

-r 

0102 

R 

0104 

*c 

0106 

eig 

0108 

A 

0110 

^cnct 

0112 

0114 

Mc 

0116 

Mig 

0118 

MP 

58 


Location 


Content 


0120 

0122 

0124 

0126 

0128 

0130 

0132 

0134 

0136 

0138 

0140 

0142 


«r 

vo 

*o 

fl 

f2 

1 

2 

0.18,  tube  coefficient  of  friction 
32.17,  g 
777.5,  J 
3 
12 
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Closed  Brake  Data 
Storage 


Location 

Content 

0070 

S01  , 

ln2 

0072 

S11  , 

In2 

0074 

S21  , 

in2 

0076 

°*1  > 

radians 

0078 

fii. 

radians 

0080 

S02  , 

in2 

0082 

S12  , 

in2 

0084 

S22  , 

In2 

0086 

<*2  , 

radians 

0088 

CM 

radians 

C090 

s03  , 

in2 

0092 

S13  , 

In2 

0094 

S23  , 

in2 

0096 

°*3  , 

radians 

0098 

7^3  , 

radians 
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Additional  Closed  Brake  Data 
Storage  for  Free-Periphery  Designs 


Location 

0052 

0054 

0056 

0058 

0060 

0062 

0064 

0066 

0068 


Content 


4 

4, 

4 

4. 


n  , 

21 


31 


12 


^22 

s  . 


4 

4. 


32 


13 


4, 


23 


33 


radians 


Open  Brake  Data 
Storage 


Location 

Content 

0090 

radians 

0092 

I  , 

inches 

0094 

Q  , 

inches 

0096 

zm  » 

inches 

0098 

*1  . 

radians 
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General  Working  Storage 


Location 

Content 

0400 

ft  -  v0 

aav  init 

0402 

init 

0404 

A  pb  init 

0406 

A  8  pb  init 

0408 

aav  init 

0410 

ab  init 

0412 

ab  init  “t 

0414 

Mt  ab  init  “  A  pb 

K  a 

0416 

A  g  pb  init 

ab  init 

0418 

^r"  A  *  Pb  init 

ab  init 

0420 

P2 

0422 

open 

0424 

open 

0426 

4*  av 

0428 

n^av 

0430 

^  av 

0432 

Tr)2 

0434 

< 

o 

to 

2gj 
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Location 


Content 


0436 

0438 


M 


Mt 

1~ 


1  +  v. 


0440 


0442 


+  MP 


3 


1  +  Ii) 

vo  / 


7 


0444 

<T  -  1)  E, 

0446 

Us  R 

0448 

EgaT 

0450 

12 

0452 

t/  Q2  -  X2 

7 T 

0454 

1  -  COS 

0520 

(Po/Poo  )2 

0522 

(P/Pc»)3 

J 
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Closed  Brake  Working  Storage 


Location 

0400 

0402 

0404 

0406 

0408 

0410 

0412 

0414 

0416 

0418 

0420 

0422 

0424 

0426 

0428 

0430 

0432 

0434 

0436 

0438 


Content 
aps  -  2  r 

bps  -  2  r 

J 

cps  «•  2  r 

J 

dps  “Sr 

J 

space 


j 

j 

j 

j 


n2  j 

n2j 

nU 


cos  OL  j 


cos  {3  j 


j- 
1  J- 


SQ 


j-1 


Sqj  +  Sij 


£ 


Ij 


S2J 


Ot 


1 

i 


€ 


2J 


nU 

n2j 


rj 

nlN  (1  ~  2  r  ) 
J-l  J 
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Location 


Content 


0440 

0442 

0444 

0446 

0448 


4 

4 

4. 


ij 

2j 


3j 
cos  4 


Ij, 


cos 


4 


Ij 


-  cos 


J 


2j 


,/2  - 


n2j  r  j  Jcos^ij-coaT/2j  +  (ff/2“'lAj)  ’ 


Open  Brake  Working  Storage 


Location 

Content 

0450 

I2 

0452 

2^  Q2  -  I2 

n 

0454 

1  -  cos  CX  . 

0456 

open 

0458 

'•mf£ 

0460 

0462 

a. 

0464 

0466 

0468 

a 

0470 

bl 

0472 

b2 

0474 

b3 

0476 

8 1 

N  | 

1 

a 

b 

rN 

0478 

to 

0480 

_  3 

0482 

(1  +'1)Pao 

0484 

_£ 

V 

0486 

1  "  jL 

V 

04HB 

66 


Location 

Content 

0490 

rN 

0492 

1  +  /lE_ 1 

1*n7 

0494 

n21  “  n(€2i> 

0496 

P»/  P« 

0498 

rl 

6  6  (16  Sept  60) 


P3  ~  CARRIAGE  SETUR* 

=  CORDIT I  OH  A  L  STOP  CODE 


PREPARED  FOR:; 


LGP-30  COOING  SHEET  Dictator  Format 


jJOS  NO.-  PROGRAM  NO.  PAGENO. 


SUPERSEDED 
PROGRAM  NO. 


II 

PROGRAM  CHECKED  BY:- 

DATE  PREPARED:. 

DATE  RUN:  ! 

L 0 C  A  T  I  ON 

(M*  «■  *0*  ClkUMMt) 


00  Lo  25  50  75 


01  26  51  76 


02  27  52  77 


63  78 


OH  2  3 


05  30  55  80 


06  31  56  81 


07  32  57  82 


08  33  58  83 


09  34  59  84 


10  35  60  85 


II  36 


12 


86 


62  87 


0104 


3  0106 


0000 


0122 


0000 


0000 


0000 


0000 


0114 


0000 


0000 


0438 


0440 


0000 


0000 


0000 


3  0440 


0000 


0438 

0000 


0  R  DON  For*  ||65  (16  Sent  60) 


38  63  88 


14  39  64  89 


40  65  90 


66 


42  67  92 


I 

oo  i  9  <J 


19  44  69  94 


II 

« 


70  95 


96 


72 

97 

73 

98 

79 

99 

I 


1391  0272 


1392 


1393  0336 


1394  0368 


13^5  0650 


0114  0518 


0116  0000 


0518  0518 


0122  0000 


0132  OOOO 


0136  0000 


0138  0434 


0118  0520 


0116  0436 


0132  0438 


0118  0440 


0120  0000 


0000  0440 


0130  0438 


0000  0000 


0140  0442 


0440  0432 


IS* 

m 


ISC 


_ 


MUL'/WAl 


ffffiiwas£g£52agl 

_  j 


0440  0000 


0442  0442 


=  CARRIAGE  0£Tij»N 
=  CONDITIONAL  STOP  CODE 


107U  99 


0  R  0  0  W  form  1165  (  !  6  Sept  60) 


S3  =  C  A  R  R  I  A  Gf  Rl  TURN 
■  COHO  I  T  I  ON  A  l  $TOP  CODE 


PREPAREO  FOR 


LOP-JO  CODING  SHEET  (Dictator  Format 


JOB  NO.  [PROGRAM  NO. 


PROGRAM  PREPARED  BY:  IPROGRAM  CHECKED  BY:  IDA TE  PREPARED:  IDATE  RUN: 


£S3  =  CARRIAGE  RETURN 
’  =  CONDITIONAL  SIOP  CODE 


ORDOW  Form  1166  (16  Sept  60) 


?  3  4  8  73  38 


1  12  4  49 


8009  0024 


0R00W  'orm  1163  (16  Sept  60) 


CARRIAGE  3E  T  U  H  H 
•  ^  CONO 11  I C N A L  STOP  COUP 


CE3  =  CARRIAGE  RETURN 
•  =  CONDITIONAL  STOP  CODE 


CED  =  CARRIAGE  RETURN 
•  -  CONDI T  I  ON  A  E  S  TO  ^  CODE 


£S3  =  carriage  return 

■  r-  CONDI  T  i  ON  A  l  STOP  CODE 


OROOW  for*  1155  (16  Sept  60) 


E3  "  C  *  R  ft  I  A  GF  Rc  TURN 


0  R  DO  W  Tor-  1  1  6  6  (16  Sect  60) 


gg  z  CARR  t  A Gf  Rf  TORN 

CONDITIONAL  S  '  0  R  CODE 


ORDOW  Form  )  i  6  5  (1C  S  e  t>  t  60), 


E3  =,  C  A  R  R  f  A  OF  Kf  TU“N 


0  R  0  0  W  F0"S  it  6  b  {  I;  6  S  e  P  5 


E3  =  CARR  I AGt  Rf  TURN 
•  CONDI  T  I  ON  A;  S  TO  !'  C  0  0  F 


0  R  D  0  *  Form  1166  (It  Sept  60! 


1553  -  CARRIAGE  RE  TURN 
v  -  CONDI T I  0  N  *  L  S  T0P  lOOt 


0  R  DO  w  form  1166  (16  Sept  60) 


gg  Z  CARRIAGE  RETURN 
■  =  CONDITIONAL  STOP  CODE 


[S3  --  CtRRU&l  Rf  TURN 

CORD  I  [  IO"M  ri  toe  mill. 


HlillUtonw 


PREPARED  FOR:-  JOB  HO, 


PROGRAM  PREPARED  BY:  PROGRAM  CHECKED  8Y:  DATE  PREPARED:-  DATE  RUN: 


problem:. 


IP RO G R AM  NO.  j 

IPAGE  NO.  1 

HH 

1400  25  50  75 


01  26  51  76 


14  02  27  52  77 


03  28  53  78 


OR  29  54  79 


05  30  55  80 


56  81 


07  32  57  82 


53  83 


0534 


3  0536 


09  34  59  84 


35  60  85 


36  61  86 


12  37  62  87 


13  38  63  88 


14  39  64  89 


40  65  90 


66  91 


1417  42  67  92 


68  93 


19  44  69  94 


95 


96 


72  I  9^ 


20 


1421 


22 


23  48  73  96 


DD 


74  99 


0510  0546 


0504  0548 


8008  0004 


8009  0C16 


8009  0006 


8009  0003 


8009  0003 


18  '  0500  0508 


8009  0016 


8009  0010 


8009  0003 


8009  0003 


0510  0518 


8009  0016 


8009  0014 


8009  0003 


8009  0003 


0000  1418 


0546  0552 


0536  0536 


0000  1425 


0546  0548 


8009  0024 


8009  0024 


18  '  0536  0536 


0  R  DON  for*  ||65  M6  S»Pt  60) 


£S3  =  CARRIAGE  RETURN 
■  =  CONDITIONAL  STOP  CODE 


84 


PREPARED  FOR: 


LGP-30  CODING  SHEET  (Dictator  Eormat 


JOB  NO.,  [PROGRAM  NO. 


PROGRAM  PREPARED  BY:-  PROGRAM  CHECKED  BY:-  DATE  PREPARED:-  DATE  RUN: 


MMI 


OR  DOW  Foil"  lies  i  I  6  Sept  60) 


tS3  =  CARR  I  AGE  R£  TURN 
-  -  CONDI T 1  ON  A l  S  TOP  CODf 


OR  DOW  For*  I  I  6  S  (  i  6  S  •  P  t  60) 


gg  =  CARRIAGE  RE  TURK 
'  =  CORO  I  T  i  ORAL  STOP  CODE 
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OR0OW  For*  M65  ( | 6  Sept  60) 


£3  =  CARRIAGE  REtURR 
'  -  CORD  I T ! CR A l  STOP  CODE 


PP r  PA  R E  D  FOR: 


PREPARED  FOR 


JOB  HO.  PROGRAM  HO.;  PAGE  NO. 

aa 


PROGRAM  CHECKED  BY  :■  DATE  PREPARED:-  DATE  RUN:  SUPE 

PROG 


PROBLEM: 


REMARKS 


00 

1525 

50 

75 

01 

26 

SI 

76 

02 

27 

52 

77 

03 

28 

53 

78 

04 

29 

mi 

79 

05 

30 

55 

80 

06 

31 

56 

81 

OOOC 


0000 


0000 


2  0004 


0008 


000 


SSHBh 


09  39  59  84 


60  8b 


66 


62  87 


38  63  68 


39  64  89 


16  40  66  90 


66 


42  67  92 


68  93 


44  69  94 


70  95 


71  96 


22  47  72  97 

23  48  73  98 


0000 


0420 


0000 


3  0430 


1  0000 


0424  0000 


0422  0430 


0740  0682 


0002  0010 


0006  0000 


549  74  99  . 


3  0000 


0002  0010 


0006  0000 


0010  0010 


0002  0000 


010  0000 


0006  0432 


0416  0000 


0410  0436 


0402  0402 


0416  0000 


04 tO  0410 


0432 


0406  JJ406 

I  * 

15  '  0426  OOOi) 
0136  nooo 


m 


F^ - 


0  ft  uO  4  ! 

o  r  »  | 

165  (16  5«Pt  50) 

80 

£53  -  CARR  ■  AGF  HE  T 0 « H 

=  CONDI T IONAI  $  TOP  f  0  Of 

0  R  D  0  w  for*  l  l J  5  (l«  Sept  60) 


Q3  =  CARRIAGE  RETURN 
>  =  CONS  I T I  ORAL  STOP  CODE 


p  ti  0  0  n  >  o'*  *16  6  (  I  6  S  n  P  I  60! 


gg|  z  C  i  S  ii  I  A  C !  F ,  T  0  <  N 

-  CONDITIONAL  STOP  CO  I) 


6b  (16  Sept  60) 


GE3  =  CARRIAGE  RETURN 

=  CONDITIONAL  STOP  CODE 


0  *  DO  W  for*  1165  (16  S  «  M  60) 


C33  =  C A  R  R  I  A  Q£  Rf  TURR 
•  =  CORD  (  T I  ORA  l  STOP  COD! 


SAMPLE  BRAKE  ANALYSES 


I.  An  example  of  a  type  1  brake  which  has  been  analysed 
by  use  of  DIPARDIP  is  the  slngle-baff le,  closed,  asym¬ 
metric  brake  pictured  in  figure  0.1  of  the  INTRODUCTION. 
When  used  with  the  standard  105mm  howitzer  at  zone  7, 
tie  weapon-break  system  produces  the  following  data. 


Location 

Number 

Parameter 

0100 

1.2593 

■% 

0102 

70.036 

R 

0104 

1350 

ec 

0106 

1231 

elg 

0108 

13.717 

A 

0110 

10.89b 

Acnct 

0112 

1431.17 

'V r 

0114 

2.83 

*c 

0116 

. 042857 

“ig 

0118 

33 

“P 

0120 

1492 

*r 

0122 

1550 

V 

0 

0070 

15.6228 

S01 

0072 

40 

S11 

0074 

53.19 

S21 

0076 

1.5707  (90°) 

<X  j 

0078 

0.6109  (35°) 

fix 

The  printed  results  take  the  following  form: 
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The  computer  print  out  for  the  following  sample  brake 
analyses  is  in  the  "+50"  system  of  notation.  The  "+50" 
system  is  analogous  to  the  powers-of-ten  notation  used  in 
scientific  work,  i .  e. , 

.0004172  =*  4.172  x  10“4 

8190000.0  «  8.19  x  IG6 

-27.05  -  -2.705  x  101 

In  the  "tSO"  system  these  numbers  would  be  written: 

4.172  x  10‘4  -  41720000 ’ +46 ' 

8.19  x  106  *>  81900000 1  -»-56  1 

-27.05  -  27050000" -51 ’ 

Notice  that  to  obtain  the  correct  "characteristic"  or 
exponent  of  the  "+50"  system,  the  usual  exponent  of  10  is 
algebraically  added  to  the  "base"  50.  The  sign  of  the 
number  always  preceeds  the  "characteristic."  The  "mantissa" 
is  always  eight  digits  long.  If  there  aren't  enough 
significant  digits  to  fill  eight  spaces,  fill  out  with 
zeros  as  was  done  in  the  above  examples. 


95 


ENERGIES 

tot 

prcj 

gas,  ktn 
rec  mass 
engrv 
gas,  them 

av  T 
av  p 

1  t ,  sec 

2  B  era! 

3  I  t,  lb  f  eec 

1  39982950*46 

2  20305754*51 

3  72108659*50 

1  79965901+46 

2  30972881+51 

3  13785520*51 

1  11994885+47 

2  55934945*51 

3  19785354*51 

1  15993180+47 

2  7142718*51 

3  25265284+51 

1  19991474+47 

2  85  58884!) +51 
5  30274556+51 

1  27980065+47 

2  11040747*52 

3  39053423*51 

1  39582952+47 

2  14039437*52 

3  49660414+51 

1  51977835+47 

2  16362256+52 

3  57876707+51 

1  63972720*47 

2  18173357+52 

3  64282947*51 

1  79965901+47 

2  19997605*52 

3  70735687+51 

1  99957376+47 

2  21605748+52 

3  76424030+51 

1  12394714+48 

2  22805603+52 

3  80951143+51 

1  14793691+40 

2  23716519*52 

3  83890268+51 
1.  17992327+48 

2  24405473*52 

3  86327241*51 

1  24389600+48 

2  25063028+52 

3  08653150+51 

1  29987213+48 

2  25296588*52 

3  89479300*51 
3828+996+49 


Btu 

38752570+53 

15848781+53 

44955162+51 

38098992+51 

22665848+52 

19786662+53 


SINGLE  BAFFLE,  ASYMMETRICAL  BRAKE 
U-BU 


19826223+53 

33447631*53 

p,  paia 
P,  lb  f  sen 
F  «,  lb  f 


B,  lb  f 
M,  lb  ra 
X  y,  lb  i'  sec 


dP/dt,  lb  f 
K,  Btu 
F  y,  lb  f 


dM/dt,  lb  n/aec 
dH/dt,  Btu/aee 
P- index 


35474565+53 

I3879082+5I 

31448244*54 

32359190+53 

26533619+51 

28686461+54 

29542047+53 

38081775+51 

26l89769*54 

26994444+53 

48629245*51 

23930607*54 

24686356+53 

58270820+51 

21884484+54 

20695579+53 

75167907+51 

10346655*54 

15979149+53 

95583666*51 

14165535+54 

12420617+53 

11139794+52 

11010892*54 

97162547+52 

12372832+52 

86154715+53 

70692786+52 

13614820*52 

62669238+53 

48177814-52 

14709683+52 

42709686+53 

30997101+52 

15561036*52 

27478967+55 

20328933*52 

16146742*52 

10021623*53 

11901330*52 

16615798+52 

10550543+53 

44238465*51 

17065477+52 

39217469*52 

20072478+51 

17222490+52 

17794276*52 

64627921+49 


4066C46O+54 

20408599+49 

11609319+51 

44387102+54 

3919455+49 

22194354*51 

40523923*54 

56499941+49 

31853945+51 

57028270*54 

72454126+49 

40676496+51 

33062274+54 

87173520+49 

48741304+51 

20508126+54 

11332048+50 

62875071+51 

21918599*54 

14561029+50 

79952070+51 

17037361+54 

17130521+50 

93180103+51 

13527786*54 

19185249+50 

10349399+52 

96969294+55 

21315186+50 

11388274*52 

66005506+53 

23261206+50 

12504005+52 

42518724*53 

24876549*50 

13032959+52 

27885197+53 

25973036+50 

15506131+52 

16325055+53 

26927360+50 

13898477*52 

60681950+52 

27907714*50 

14272944*52 

27533417*52 

28290306+50 

14405952+52 

83646164+49 


33129140+54 

17535141*52 

27711255*54 

30219742*54 

33371841+52 

25277654+54 

27589603+54 

47688445+52 

23077644+54 

25209690+54 

60642997+52 

21006938+54 

23054202*54 

72375959*52 

19283955+54 

19327278+54 

92662668+52 

i6i£6526+54 

14922677+54 

11663817+53 

12482247+54 

11599421+54 

13472762+53 

97024708*53 

90738584+53 

14847499+53 

75899543*53 

66018889*53 

16193063*53 

55222268+53 

44992509+53 

17337590+53 


28947709+53 

18210955+53 

24213647*53 

18904872+53 

1&753261+53 

1580117*53 

11114466*53 

19101670*53 

92968247+52 

41313651+52 

19561780*53 

34557268*52 

18745374*52 

19604626*53 

15679786*52 
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48951311+52 
41653713*55 
64627921+49 
45076975+52 
37637816+55 
64627921+4? 
41541370+52 
34041450+55 
64627921+49 
38312137+52 
30617422+55 
64627921+49 
35  36 0255  *52 
27924312+55 
64627921+49 
30187012+52 
22988940*55 
64627921+49 
23936487+52 
17203465+55 
64627921+49 
19094759+52 

13090496+55 

64627921+49 

15319668+52 

99846388+54 

64627921+49 

11517174+52 

70305320+54 

64627921+49 

81651219+51 

46059104+54 

64627921+49 

54973706+51 

28518557+54 

64627921+49 

2RHBS 

64627921+49 

23293249+51 

98524659+53 

64627921+49 

95870265+50 

33075056+53 

64627921+49 

47186572+50 

13834632+53 

64627921+49 


II.  A  hypothetical  triple-baffle,  closed,  symmetric  brake, 
which  is  used  with  the  105ram  howitzer  is  analysed  in  this 
example.  Figure  0.2  of  the  INTRODUCTION  illustrates  this 
type.  Since  the  thermodynamic  data  for  this  example  is 
identical  to  that  of  example  I,  the  contents  of  locations 
0100-0122  are  not  written  here,  only  the  brake  data  being 
given. 


Location 

Number 

Parameter 

0070 

16 

S01 

0072 

50 

Sll 

0074 

100 

S21 

007b 

2.3562 

(135°)  °<1 

0078 

1 . 5707 

(90°) 

0080 

20 

s02 

0082 

50 

S12 

0084 

100 

S22 

0086 

2.3562 

Ot  2 

0088 

1.5707 

^2 

0090 

25 

s03 

0092 

50 

S13 

0094 

100 

S23 

0006 

2.3562 

Ot  3 

0098 

1.5707 

ft  3 

The  computer 

print  out  for 

example  II  follows 
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ENERGIES 

tot 

proj 

gas,  kin 
rec  mass 
engrv 
gas,  therm 


Btu 

38732570+53 

158^8701+53 

44955162+51 

36098992+51 

22665848+52 

19786662+53 


av  *  19826223+53 

av  P  33447631+53 


1 

2 


t»  sec  P,  poia 

®  cml  P,  lb  f  sec 

I  z,  lb  f  sec  P  a,  lb  f 


2  20397554+51 

3  12553620-51 

1  79925189+46 

2  38995440+51 

3  23999640-51 

1  11988778+47 

2  55967319+51 

3  34444937-51 

1  15985037+47 

2  71468532+51 

3  43985117.51 

1  19901296+47 

2  85638383+51 

3  52705915-51 

1  27973815+47 

2  11047138+52 

3  67989318-51 

1  39962595+47 

2  1+047563+52 

3  86455360-51 

1  51951372+47 

2  16371726+52 

3  10075936-52 

1  63940151+47 

2  16183875+52 

3  11191218-52 

1  79925188+47 

2  20009170+52 

3  12314595-52 

1  99906486+47 

2  £1610253+52 

3  13304896-52 

1  12388404+48 

2  22898849+52 

3  14093036-52 

1  14786160+40 

2  23730246+52 

3  14604717-52 

1  17983167+48 

2  24419599+52 

3  15028977-52 

1  24377182+48 

2  25077534+52 

3  15433902-52 

1  29971945+48 

2  25311229+52 

3  15577729-52 
99445308-49 


39513170+53 

13887115+51 

76693976+54 

32394413+53 

26548977+51 

71783076+54 

29575004+53 

38103817+51 

65535522+54 

27023828+53 

48657391+51 

59882348+54 

24713227+53 

58304548+51 

54762266+54 

20718106+53 

75211414+51 

45909440+54 

15996542+53 

95638990+51 

35446885+54 

12434137+53 

11146241+52 

27552919+54 

97266307+52 

12379993+52 

21553773+54 

70769733+52 

13622700+52 

15681930+54 

40230254+52 

14710197+52 

10687386+54 

31030841+52 

15590054+52 

68761526+53 

20351060+52 
I6156069+52 
4509610c ■ 13 

11914264+52 

16625415+52 

26400972+53 

44286638+51 

17073354+52 

98135167+52 

20094327+51 

17232450+52 

44527202+52 

16154474+50 


105  !♦(  HOWITZER  AT  ZONE  7 
WITH  A  TRIPLE -BAFFLE,  CLOSED, 
SLWETRIC  BRAKE 


B,  lb  t 
M,  lb  m 
1  y,  lb  f  sec 


dP/dt,  lb  t 
H,  Btu 
p  y#  lb  f 


dM/dt,  lb  b/mo 
dH/dt,  Btu/aec 
P -index 


48713426+54 

20408599+49 


33165201+54 

17535141+52 


44435416+54 

39194355+49 


30252636+54 

33371841+52 


40568032+54 

56499941+49 


27619633+54 

47688445+52 


37068584+54 

72454126+49 


25237131+54 

60642997+52 


33899134+54 

87173520+49 


23079296+54 

72375939+52 


28419026+54 

11332048+50 


19348315+54 

92662668+52 


21942457+54 

14561929+50 


14938920+54 

11663817+53 


17055906+54  11612046+54 

17130521+50  13472762+53 


13342293+54 

19185249+50 


90837350+53 

14847499+53 


97074843+53  66090750+53 

21315186+50  16193063+53 

66157440+53  45041482+53 

23261206+50  17337590+53 


42565004+53 

24876549+50 


28979218+53 

16210955+53 


27915549+53 

25973036+50 


19005537+53 

18753261+53 


16342824+53 

26927360+50 


11126564+53 

19101670+53 


60747981+52 

27907714+50 


41358600+52 

19561780+53 


27563388+52 

26290306+50 


18765778+52 

19684626+51 


OOOOOOOO+49 


48976245+58 

41674931+55 

16154474+50 

45099936+52 

37656988+55 

16154474+50 

41562530+52 

34058769+55 

16154474+50 

38331653+52 

30833120+55 

16154474+50 

35378260+52 

27938537+55 

16154474+50 

30202368+52 

23000650+55 

16154474+50 

23948680+52 

17292269+55 

16154474+50 

19104486+52 

13097165+55 

I6154474+50 

15327471+52 

99897247+54 

I6154474+50 

11523041+52 

70341132+54 

I6154474+50 

81692811+51 

46c02>65+54 

16154474+50 

55001708+51 

28332982+54 

16154474+50 

37673115+51 

17791965+54 

I6154474+50 

23305114+51 

98574846+53 

16114474+50 

95919099+50 

33091903+53 

16154474+50 

47210609+50 

13841679+53 

I6I54474+10 
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III.  The  single-baffle  open  brake,  pictured  in  figure 
0.3  of  the  INTRODUCTION,  has  been  analysed  using  DIPARD1P. 
When  this  brake  is  used  with  the  XM103  cannon  at  zone  11, 
the  following  system  data  are  used  in  the  analysis. 


Location 

Number 

Parameter 

0090 

1.4388  (82. 5n) 

0092 

2.21 

I 

0094 

7.86 

Q 

0096 

5.85 

z» 

0098 

1.9635  (112.5° 

)  oi  x 

0100 

1.2402 

7 

0102 

67 

R 

0104 

1727.5 

ec 

0106 

1231 

eig 

0108 

13.717 

A 

0110 

10.8984 

Acnct 

0112 

17  lu 

T/t 

0114 

4.46 

«c 

0116 

0.042t  '.7 

Mig 

0118 

28.5 

MP 

0120 

1430 

Mr 

0122 

2200 

Vo 

The  results  of  the  analysis  of  this  example  follow. 
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ENERGIES 

tot 

proj 

gas,  kin 
rec  mass 
engrv 

gas,  therm 

av  T 
av  p 

1  t,  aec 

2  B  cml 

3  I  if  lu  T  sec 

1  586124994*6 

2  56299054+51 
5  12910970-50 

1  77224996+46 

2  6  95951*  19+51 
5  20125245-50 

1  13565749+47 

2  99598199+51 
5  l'5265971-50 

1  154U4999+47 

2  12718581+52 
5  6726fr762-50 

1  19506249+47 

2  15240016+52 
5  90214555-50 

1  27020748+47 

2  19659241+52 
5  15&769&7-51 

1  58612499+47 

2  24990755+52 
5  21440072-51 

1  50196246+47 

2  29154762+52 
5  29107960-51 

1  61779996+47 

2  52559622+52 
5  57207181-51 

1  77224096+47 

2  55607309+52 
5  45242451-51 

1  96551245+47 

2  50471564+52 
5  50511224-51 

1  11969074+48 

2  40750204+52 
5  52321012-51 

1  14286624+48 

2  42229020+52 
5  52350629-51 

1  17375624+48 

2  45456576+52 
5  51162351-51 

1  255  5  5624+48 

2  44627425+52 

3  48521706-51 

1  28959575+48 

2  45045505+52 

3  47  357602-51 
1105  1  579+50 


Btu 

77574071+55 

27574571+53 

14217041+52 

63701100+51 

39435274+52 

43996456+53 


XM103  CANNON  AT  ZONE  U 
WITH  SINOUE-BAFFIE, 
GERMAN,  OPEN  BRAKE 


27235043+53 

57660201+53 


p,  psia 
P,  lb  f  sea 
F  s,  lb  f 

65408216+53 
24713203+51 
92911659+54 
59664069+53 
47245970+5  i 
85602916+54 
54471270+53 
67800710+51 
78959169+54 

49772518+53 

86589616+51 

73718547+54 

45516853+53 

10375748+52 

60116015+54 

38158653+55 

13384456+52 

58093190+54 

29462451+53 

17019702+52 

46158612+54 

22901209+55 

19835604+52 

37297507+54 

17914888+53 

22031162+52 

30745566+54 
13034578+53 
24242656+52 
22258020+54 
38850542+52 
26192176+52 
14541798+54 
57152642+52 
27743715+52 
85310077+55 
37402607+52 
20751017+52 
51517815  '53 
21943744+52 
29586221+52 
27184625+53 
01567184+51 
3O38336I+52 
86651005+52 
37009756+51 
30666501+52 
36556079+52 

00000000+49 


B,  lb  f 

M,  lb  m 

I  y,  lb  f  sec 

dP/dt,  lb  f 
R,  Btu 

F  y,  lb  f 

89720448+54 

31988020+49 

61083708+54 

38990108+52 

81841205+54 

61432428+49 

55719341+54 

74203662+52 

74718253+54 

88556846+49 

50869874+54 

10603722+53 

68272964+54 

11556310+50 

46401779+54 

13484220+53 

62435467+54 

13663398+50 

42507477+54 

16093087+53 

52542195+54 

17761619+50 

35635751+54 

20603925+53 

4o4l3644+54 

22824 06  8+50 

27514522+54 

25934978+53 

31413589+54 

26850028+50 

21387082+54 

29957240+53 

24573855+54 

30070565+50 

16730434+54 

35014025+53 

17879256+54 

33408984+50 

12172602+54 

36005942+53 

12164085+54 

36459137+50 

02957*+58+53 

38550845+53 

78396280+53 

38990992+50 

53373959+53 

40492808+53 

51414895+53 

40709602+50 

35004422+53 

41698648+53 

30100234+53 

422*^591+50 

20492919+53 

*'2651233+53 

11168570+53 

43741977+50 

76174325+52 

43496423+53 

50766283+52 

44341646+50 

54562855+52 

43769575+53 

100 

dM/dt,  lb  m/sec 
dH/dt,  Btu/sec 
P- index 


79448446+52 

95906011+^3 

10355683+50 

73160362+52 

86659571+55 

10459635+5O 

67422042+52 

78379033+55 

10567587+50 

62180967+52 

70955805+55 

10797619+50 


57390035+52 
64294616+55 
10909955+50 
40993011+52 
52931116+55 
IIO98729+50 
38849140+52 
39794489+55 
11416592+50 
50990977+52 
30140346 +55 
11&73073+50 
24863968+52 
2?989232+55 
12511495+50 
10692405+52 
16187519+55 
12449074+50 
13252071+52 

10604924+55 
11770154+50 
89222852+51 
65202346+54 
10081904 +?0 

61112697+51 

40944433+54 

10020018+50 

57605166+51 

22684909+54 

90313670+49 

15559036+51 

76153995+53 

77446717+49 

76584261+50 

31853609+53 

72008581+49 


IV.  As  a  final  example  of  the  use  of  DIPAP.DIP,  a  double¬ 
baffle,  closed,  free-periphery  brake  was  chosen.  This 
is  pictured  in  figure  0.4  of  the  INTRODUCTION.  The  brake 
was  designed  to  be  used  with  an  improved  155mm  howitzer 
for  which  the  system  data  are  given  below. 


Location 

Number 

Parameter 

0052 

0.78530 

^11 

0054 

1.44835 

nJ-  21 

0056 

0.314136 

^31 

oor.a 

0.78530 

-^12 

0060 

1.44835 

'^22 

0062 

0.314136 

^32 

0070 

31.371 

S01 

0072 

139.0 

S11 

0074 

16.36 

S21 

0076 

2.3562  (135° 

)  <*1 

0078 

1.5707  (not 

used)  fil 

0080 

32.170 

S02 

0082 

139.0 

s12 

0084 

222.5 

s22 

0086 

2.3562 

CL  2 

0100 

1.2402 

0102 

67 

R 

0104 

1727.5 

®c 

0106 

1500 

eig 

0108 

29.823 

A 

101 


Location 

Number 

Parameter 

0110 

19.30 

Acnct 

0112 

6707.1 

¥r 

0114 

24.08 

Mc 

0116 

2.714  10“3 

“ig 

0118 

95 

“p 

0120 

4900 

Mr 

0122 

2700 

Vo 

The  print  out  of  this  analysis  follows. 
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ENEKJIBo 
to  l 
proj 

gas,  kin 
rec  mass 
engrv 

gas,  them 

av  T 
av  p 

1  t,  sec 

2  B  cml 

3  I  2,  it  *  SCC 

1  62322650+46 

2  190611543+52 

3  11966464-52 

1  12564529+47 

2  37976446+52 

3  22877171-52 

1  18646794+47 

2  54504832+52 

3  32833939 -52 

1  25129060+47 

2  69600981+52 

3  41927923-52 

1  31411324+47 

2  83400559+52 

3  50240845-52 

1  43975853+47 

2  10758^64+53 

3  64809440-52 

1  62822651+47 

2  13680485+53 

3  82411819-52 

1  81669443+47 

2  15943916+53 

3  96046818-52 

1  10051623+48 

2  17708711+53 

3  10667602-53 

1  12564529+48 

2  19486316+53 

3  11738636-53 

1  15705662+48 

2  21053345+53 

3  12682623-53 

1  19475021+48 

2  22300477+53 
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APPENDIX  I 


Momentum  and  Energy  Distribution 
at  Start  of  Gas  Discharge  Period 


In  a  frame  of  reference  which  is  stationary  with 
respect  to  the  ground,  the  assumed  velocity  distribution 
m  the  gas  behind  the  projectile  within  a  gun  tube  of 
uniform  cross  section  and  length  x  is  shown  below  in 
figure  1.1 

v 

Figure  1.1 


11 


x  •  £  *  ~r 

Let  x  -  n  ,  then 

IT  1 


v 


/v.  +  vr 


) 


X 


-  V 


r 


Now,  momentum  conservation  requires  that 

1  2a  Mr  vr  -  PT  +  Mp  v0  or 

8  8~ 

b  vr  "  pr  +  Mp  v„/g 

Mr''* 


However,  PT  ,  the  momentum  taken  by  the  charge  and 
igniter  gases,  is  found  by  averaging  the  momentum  dis¬ 
tribution,  i  ■  e . , 
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1  ,3a 


'  T 


«T 

gL- 


v  dx  or 


Vo 


Mtv0 


'*  UVr) 


1  +  V 


r 

~ 


h 


V,. 


«Tv, 


1  -  vr  ' 


1  4 


Substitution  of  1.3b  into  1.2b  gives 

Mt  +  Mp 
vr  2~ 


My  v  Mr 
2~ 


Since  the  projectile  leaves  the  muzzle  with  velocity 
it  carries  off  an  energy  given  by 


KEP 


1  Mp  v- 

2 

g 


In  engineering  units  the  kinetic  energy  of  the  projectile 
in  Btu  is 


I  5 


KEP  -  Mpv‘ 


2gJ 


Similarly,  the  kinetic  energy  of  the  recoiling  parts  is 

,2 


1  6 


KER  -  Krv‘ 


with  vr  being  given  by  equation  '  4 


2gJ 


The  kinetic  energy  of  the  gas  itself  may  be  found  by  Inte¬ 
gration  of  an  elemental  mass  of  gas  o-’er  the  entire  internal 
volume  of  the  gun  tube.  Thus 


1  .  7 


KEG  -  1  M  v?  /J./VJ2  dn 
2  _! —  u  \v0/  ! 


g 


Substitution  of  1.1  for  v_  and  integration  yields, 
in  Btu  v0 

1.8  KEG  -  1  Mt  v2  /l+vr\3  +/vr\2  -  vr  /l+vr\]  . 

^  \  v7/  \v7/  v7  V  v~J 

After  in  assumed  initial  shear  deformation  in  the 
ammunition  band  as  the  projectile  leaves  the  chamber, 
during  which  relatively  little  energy  is  expended,  work 
is  performed  in  forcing  the  projectile  thru  the  tube 
against  frictional  forces.  The  energy  thereby  expended, 
called  the  engraving  energy,  is  computed  as  follows. 

The  minimum  pressure  which  mRy  exist  between  the  band 
and  the  walls  of  the  tube  to  effect  a  seal  is  the  pressure 
at  the  base  of  the  projectile,  p.  Using  .18  as  the  effec 
tive  coefficient  of  friction  between  the  band  and  walls, 
the  frictional  force  is  given  by 


•A8  Acnct  P  • 

The  work  performed  by  this  force  in  Btu  is, 

f  o 

1.9  KEB  -  .18  Acnct  p  d£  ,  where  ^ 

travel  of  the  projectile  during  band  deformation. 
However,  KEP  -  A  /  p  d  ^  . 

Hence, 


then, 
is  the 


■  ’0  KEB  -  -l«Acnct  KEP 

- x — 

The  sources  of  energy  for  the  sinks  mentioned  above 
are  the  energies  obtained  by  the  combustion  of  the 
igniter  and  charge.  Thus 

l.Ha  7/ig  -  eig  Mig 

b  "  ec  Mc  >  *here  e^g  and  ec  are  the  heats 

of  explosion  for  the  igniter  and  charge,  respectively. 

The  total  energy  is,  therefore, 

i.i2  -  yle  *  . 

At  the  start  of  the  gas  discharge,  the  thermal  a lergy 
remaining  in  the  gas  (assumed  to  be  completely  combusted) 
can  be  round  by  subtracting  all  of  the  above  kl/ietic 
energies  from  the  total  energy.  Thus, 

Eg  -  ?/•-  S  KE’s 
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1.13 


This  energy  can  be  expressed  in  terms  of  the  average 
initial  temperature  of  the  gas,  Tg  «  Tav  init* 


1.14 


Eg  -  Tg  MtCv 


However,  Cv  -•*  R _  Substituting  the  latter 

('£  ~1)  j 

expression  into  equation  1.14  and  solving  for  Tg  gives 

1.15  Te  -  Eg  ('*  -1}  J 

R  Mt 

At  te  the  gas  behaves  in  a  nearly  ideal  fashion  so  that 

1.16  pV  -  RT  where 

IS 

1.17  V  -  fJr  . 


From  the  last  three  relationships,  one  can  solve  for 
the  average  initial  pressure  in  the  gas  as 


1.18  Pg  -  12  Eg  (t-1)  J  . 

*Vr 

Figure  1.2  shows  graphically  the  energy  distribution 
at  the  start  of  gas  discharge  for  the  standard  105mm 
howitzer  (zone  7).  Values  for  the  energy  of  all  sources 
and  sinks  (including  some  for  which  no  equations  are  given 
here)  have  been  computed.  Some  are  seen  to  be  negligible 
and  were  not  considered  in  subsequent  analyses. 
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PERCENT  DISTRIBUTION  OF  FNERGY 
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APPENDIX  II 


Derivation  of  the  Relationships: 
4>av  (/3>  *  av 


Having  determined  the  average  pressure,  temperature, 
and  specific  volume  in  the  gas  within  the  gun  tube  at  t„, 
it  is  desirable  to  compute  their  local  values  anywhere 
along  the  tube  length— and,  in  particular,  at  the  breech— 
as  a  function  of  the  average  values  and  the  muzzle 
velocity. 


As  before,  we  examine  a  tube  of  uniform  cross  section 
of  length  L  and  total  internal  volume  ,  and  we  define 
a  dimension  x  such  that 


2.1 


2.2 


"  £  ♦  ,  where 

A 


-  I 


Vo. 

A 


For  convenience,  let  us  also  define  some  dimension¬ 
less  variables: 


2.3a 


2.4 


<P  “  P 

pb 

-  T 

Tb 

f  -  V 

1'r- 


and 


Since  the  average  temperature  within  the  tube  was 
obtained  from  the  thermal  energy  remaining  in  the  gas,  it 
is  necessary  to  require  that  the  enthalpy  integrated  over 
the  volume  of  the  tube  equal  the  thermal  energy.  Thus, 


2.5 


2.6 


Sg  -  Mt  Cv  Tav 


-  Cv  A  /0  T  dx 


or 


Tav  -  AL 


fL  T 
0  Y~ 
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If  we  regard  T(rj  ),  V(r|  ),  then  2.6  can  be  expressed 


Tav  “  Vav  /qX  (r\)dr) 


However,  in  an  ideal  gas  p  -  RT  ; 

T2  r~ 


therefore. 


Pav  "  P<1)  d*)  •  or 


in  non-dimensional  terms 


$  av  “  ^  ( *1  )  d  ^  • 


Sutton  (reference  book  5,  p.  81)  has  shown  that 
isentropic  decay  and  momentum  conservation  within  the  tube 
demands  that 

2.9  1_  -1+  ^  M2  ,  where  M  is  the  local  mach  number. 

4> 

Also  for  isentropic  expansion: 

2.10a  l/  -  fp 


\ 

-  yi  a  . 

M  can  be  written  as 


M  «  v_  -  v 

a  ab(a/ab) 


a  *b(a7ab)  ^  ^  (If  -1)/ 2 Tf 
As  an  approximation  to  equation  1.1,  we  may  state 


By  further  defining 


vc  -  OL  ,  2.11  can  be  written  as 


a  (t  -1) 

t  “TT 
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Substitution  of  the  latter  expression  for  M  into  2.9 


gives 

2.15 


n  -/  i-  4> 

^  ot  2  §x/* 


1/2 


This  equation  expresses  implicitly  as  a  function 
of  ^  ;  and,  unfortunately,  does  not  permit  a  convenient 

explicit  formulation.  However,  numerical  evaluation  per¬ 
mits  one  to  establish  $  (  n ) .  Then,  numerical  evaluation 
of  equation  2.8b  determines  the  relationship  a  (  ot.  ). 


2.16a 

b 


For  mass  conservation, 

m  -  a  jL0  r dx 

1  “  Vav  ft  1 

aV  0  YVa  ) 


Division  of  2.16b  by  Vb 


or 

d 

gives 


2.17  1  -  /  d  h  ,  which  with  2.10b  give 

0  yTTp- 

2-ie  • 


This  expression  also  must  be  evaluated  numerically. 

Since  the  gas  behaves  nearly  ideally, 

2.19a  pV  -  RT  and 

12 

b  pbVb  -  RTb  . 

~n~ 

Thus, 

2.20a  <f>  ^  and 

b  ^av  ^av  *  ^av 

If  we  define  a  parameter  /3  such  that 
2.21  ^3  -  v0  f  it  follows  that 
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Having  determined  <j)  ay,  ^ay,  *nd  i)*av  as  functions 

of  <X  ,  one  can  with  equation  2.22  determine  these  as 
functions  of  £  . 

An  optimum  quadratic  fit  to  the  numerical  values 
O  av^/3)  and  ^  av(/S)  was  performed  using  a  least  squares 

criterion.  These  quadratics  were  used  in  DIPARDIP  to  cal¬ 
culate  av  and  g  av  to  a  degree  of  accuracy  warranted 
by  the  data.  ' 

After  projectile  ejection,  the  gas  rapidly  becomes 
sonic  at  the  muzzle  and  equation  2.9  becomes 


APPENDIX  III 


State  Variables  During  the  Discharge  Period 


In  appendix  III  are  derived  the  equations  governing 
the  dimensional  and  non-dimensional  state  variables  as 
functions  of  time.  The  relationship  n(£  )  is  also  derived 
for  one-dimensional  isentropic  flow. 


For  one-dimensional  gas  discharge  from  a  tube  or 
choked  nozzle,  the  isentropic  mass  rate  of  discharge  is 
given  by  (p.  61,  reference  book  5). 


3.1a 


b 


M 


A 


Pb  g  nr 

ab 


f+l 

2(15-0 


A  g~* 


(V») 

2(15-0 


Pb 


<*b 


The  average  specific  volume  of  the  gas  within  the 
tube  is  the  total  internal  volume  divided  by  the  mass 
present  within  the  tube.  The  mass  present  at  time  t 
equals  the  initial  total  mass  minus  the  mass  discharged. 
Thus, 

3-2  - _ Ml _ . 

*  T  ^0  ®dt 

Also, 

3‘3  V*v  -  \  av  Vb  *ad 

3,4  Vb  init  -/Pb _ \l/%. 

^b  \P'd  init  / 

From  3.3  and  3.4  one  obtains 

3,3  ^av  "  V*v  init  H  av  /  Pb  init  \ 

^  av  init  \  pb  / 

Elimination  of  Vav  by  substitution  of  the  latter 
expression  into  3.2  yields 
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3.6 


^ftv  init  K  av  /  Pb  init  \ 


^  av  init  \  Pb  / 

_  %r _ 

Mt  "  C1  /t  pb/ab  dt 
0 

As  long  as  the  gas  velocity  at  the  muzzle  remains 
sonic,  av/  av  inj  t  -  1.  (Equation  2.23  implies  that 

§  av  is  constant  for  a  choked  nozzle.) 

Therefore,  3.6  becomes 

3-7  /Pb  init^K  "  l/r/pb  init  Vav  init 

^  pb  I  ^  “  C?  7TZ  pb  Ht  . 


v  pb  1  ii _ -  c;  jj _ p^  „.  . 

Pb  init  Pb  init  ab 

Let  us  define 


f  so  that  equation  3.7  may  be 


written  as 


pb  init 


1  m  1/t  ^  init  Vav  init 

r2*  ** _ :'ct  7j 

pb  init  °  ab 


But»  Tb  init  -  / Pb 


"S  -1 


3.10 


Tb  "  Tb  init  1  ’1>/‘* 

ah  -  (  *8  g  R  Th)1/Z  • 


so  that 


ab  “  ab  init 


-1)/2K 


Substitution  of  ab  in  equation  3.12  for  that  value  in  3.9 
yields 

3.13  X1^  - 


u  _  _  n*  t  *tf+ 1 

— -  -1 -  !  $  2T  dt 

pb  init  ab  init  0 
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pb  init  Vav  init 

ITt 
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Finally, 


3.20 


<&  “  *  ( 1  -1)  /_2 

2  Tf 


Y  +1  t-/2  if 


(-br) 


T  =T) ^  +1  \  Z 


w 


To  compute  the  momentum  ra *  of  discharge  from  the 
muzzle,  one  observes  that 


P  ■  M  a! 


And  from  equations  3.1a  and  3.21  is  derived 
£  “  A  Y  ao  / _ 2_\  2("ir+-l') 

si — (t*;  y  pb  • 


From  equations  2.10c^and  2.23,  one  has 
3.22  a;  -  (l+r>  . 


gives 


Substitution  of  the  latter  expression  for  a '  in  3.22 

*b 


Y  +1 


M 


(  r  -D/2Y 


We  may  also  define  a  dimensionless  momentum  rate  of 
discharge  jjl  as 

3.24  /X-  P _ 

A  Pb  init 

Then  using  this  and  3.8,  we  have 


3.25 


-(1=1)  _l±L 


It  is  also  desirable  to  know  the  total  momentum 
discharged  to  time  t  (assuming  that  the  gas  at  the  muzzle 
remains  sonic) .  That  is 


P  ■  /  P  dt  ,  which  with  3.24  produce 
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3.27a 

P  "  A  Pb  init  /0  A  dt 

or 

b 

P  "  A  pb  init  P*  d  J 

otT 

and 

c 

P  "  ab  init  Mt  J0  /^d  J  > 

e 

from  3.19b 

Defining, 

3.28 

Mrm  /  /AdT  ,  we  can  write  3.27c  as 

0 

3.29 

p  "  ab  init  “t  4i  • 

In  a  similar  manner  to  that  in  which  M-  was  defined, 
let  us  define  a  dimensionless  mass  rate  of  discharge  V  . 


3.30a 


y  -  M  a 


b  init  , 
A  g  pb  init 


so  that 


b  M  "  A  g  Pb  init  y  . 

ab  init 

From  equations  3.1a,  3.8,  and  3,30a,  one  obtains 

3.3i  v  -  svfrW) 


V*b  i„u 


However, 


3.32 


^b _ -/£b _ 

ab  init  \Pb  init 


Y-l 

?T 


x-' 

2Y 


Therefore, 


3.33 


V  - 


K  +1  Y  +1 

zrr=Ty--rT 
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Let  us  define  the  following  constants 


Y  +i 


2TT  -I) 


3.34a 


ci  “  Y  ( 


c2  -  2  Y 


Y  -1)  /  2  \ 

-2  VynJ 


C4  > 


Y 


-T 


^(yrr) 


Y  +i 
=TT 


-  C4 


a+  <r) 


<*-i>/2y  (l+Y) 


(Y^T7/2if 


Y+i 

2('tf  -1) 

d  C.  -  y/  2  \  .  and 

4  «(' T*U 


e  C-  -  Y  +1  . 

5 

Then,  equation  3.20  becomes 
3.35  $.  -  (Cx  7  +1)"C2 


The  cumulative  breech  force  to  time  t — or  breech 
impulse — is  given  by 

,t 

3.36  B  «  /  B  dt  ,  where 

t  o 

3.37a  B  -  A  pb  or 

b  B  “  A  Pb  init  $  * 

Thus, 

r 

3.38  Bt  -  A_Pb_init  fQ  $  d  ^  . 

CXT 

Defining  a  dimensionless  breech  impulse, 

3.39  $T-/0?  d  Y  ,  we  have 


lie 


3.40 


bt  "  A  Pb  init 


$Tmay  be  found  by  substitution  of  in  3.35  into 
3.39.  One  obtains 


3.41a 


$T-  Jq  C1 **  +1  "C2  dT 


i  -  (^r  +i)’(c2'1)l  . 


v'1^2* 


Thus,  from  3.25,  3.28,  3.34c,  and  3.41b,  one  can 
write  as 


//-T  - 


i  -  ;cxr  +d‘(c2“1> 


Further,  it  is  desired  to  have  the  cumulative  mass 
discharged  to  time  t.  This  is  given  by 

3.43  M  -  J*  M  dt  . 

Substitution  of  U  in  equation  3.30b  into  3.43  gives 


-  «T  /„  *  « 


i  . 


The  latter  may  be  expressed  as 

3.45  M  -  M  T  Vt  where 

3.46  Vt  -  /  V  d  T  ,  a  dimensionless  cumulative 
mass  discharged. 

Vt  can  be  found  by  substitution  of  V  in  equation  3.33 
into  3.46  .  Using  the  constants  defined  in  3.34,  Vt  ®»y  be 
expressed  as 


VT  -  C4 


i  -  (c^r  +D 


~<C5-1) 
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The  total  energy  rate  of  discharge  or  stagnation 
enthalpy  rate  of  discharge  Is  given  by 

3.48  '  H  -  Cp  Tav  M  . 

However, 

3.49a  Cn  -  Y  R  and 

p  Ty-nv/ 

b  Tav  -  Tay  init  <|  “Y 


From  equations  3.30b,  3.48,  and  3.49, 

Tf 

3.50  H  -  YR  Tav  init  A  *  Pb  init  J  “ 

TT-tt  - vrsn - vi  1 

that  is,  with  3.33  and  3.35 

3.51  H  ■  Y  R  ^av  init  A  %  Pb  init  * 

(*y  -iyj  *b 

.  c4  ^T+l  "(C5+2)  . 

or,  from  1.15,  3.19b,  and  3.51,  we  have 


H  ■  E, 


<*t  t  c4  cxr+i  "(C5+Z)  s 


Defining  a  dimensionless  stagnation  enthalpy  rate  of 
discharge  rj  as 

r  _  1  -<C5+2) 


y\  -  *c4  fciTT  +i 

®  "  Eg  <=*T  ^  * 


»  one  obtains  for  H 


The  cumulative  stagnation  enthalpy  discharged  is,  then, 

H  -  /*  fi  dt  or  from  the  proceeding  two  relations 


-  E,  /  n 
Z  Jo  I 


d  r  . 
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With  the  definition 

or 

3.57  rjT-  /  rj  d  7  ,  equation  3.56  becomes 

3.58  H  -  Eg  r|T 

One  can  obtain  a  computational  expression  for  n 
from  3.53  and  3.57  as  follows  ' 


3.59 


n_-  ^  c4 

,r 


[l  -  (CiT +l)“(C5+1) 


All  of  the  non-dimensional  parameters  were  evaluated 
as  functions  of  T  for  ^  -  1.26  in  the  equations: 

3.20  or  3.35  for  <£ 

3.25  for  /<- 

3.33  for  V 

3.53  for  rj 

3.41b  for  $.r 

3.42  for  ^tT 

3.47  for  v»T 

3.59  for  ^ 

Auxiliary  constants  are  defined  in  3.34a  thru  e. 
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Derivation  of  the  Relationship:  n(fe) 


In  appendicies  IV  and  V  a  relationship  is  required 
which  exists  between  the  velocity  to  throat  velocity  ratio 
and  the  area  to  throat  area  ratio  in  one-dimensional, 
isentropic  flows.  The  following  is  a  derivation  of  this 
relationship. 

In  any  of  the  books  listed  in  the  BIBLIOGRAPHY,  are 
found  one-dimensional,  isentropic  flow  formulas  involving 
the  variables:  U,  reach  number; 

£  ,  expansion  ratio; 

T 

ir-f  temperature  ratio;  and 
b 

*  ratio  of  speed  of  sound  to  stagnation 

ab 

speed  of  sound. 

The  subscript  "b"  is  used  here  to  denote  a  stagnation 
value  or  the  value  of  the  variable  at  the  breech  of  a 
tube.  The  subscript  *Vf  denotes  the  value  of  the  variable 
at  the  muzzle  or  throat.  Variables  without  subscripts  are 
evaluated  at  a  station  downstream  of  the  throat. 

The  formulas  used  in  the  derivation  are: 


€  -  s 


,  by  definition. 


3.65 


n  -  v_  ,  by  definition. 
»•  -  Tf  +1 

2  1+T*  -1  *-=T 

1  -  M2  ~T~ 

7*  I+TTHP 


from  isentropic  flow  theory. 

“  v^  ,  by  definition  . 

8  R  T 

2 

Tb  -•  1  ♦  ^  -1  M  ,  from  isentropic  flow  theory. 


T  LJ 

a*  =  (1+**)  ** 

*b 


,  from  isentropic  flow  theory. 
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1/2 

3.66  *tj  ■  ( Y  g  R  Tjj)  s  by  regarding  sound 
propagation  as  a  strictly  adiabatic  process.  Substitu¬ 
tion  of  T  from  3.63  into  3.64  gives 

3.67  Th(M2Y  g  R)  -  v2(l+  Y  -1  M2)  ,  which  by  3.66 

“2 

becomes 

3.68  M2ab2-  v2(l+|-l  M2) 

2 

Elimination  of  a^  between  3.65  and  3.38  yields 


3.69 


„  til 

n2  -/M2  “  M2  (14  Y)  * 

\  *o/  1+7  -1  M2 


3.70 


By  solving  this  expression  for  M2,  we  get 
M2  -  n2 


'  V-l  /YL-iv,2 
(1+Y  )— 


-Nr‘ 

Substitution  of  this  value  of  M2  into  3.62  produces 


3.71 


1  -  n2 


a  -  bn 


1+1 

2  +  (V-P52 

a  -  bn2 


where 


Y  -1 

a  -  (1*  Y)"T" 
b  «  if  -1 

c  -  Y  +1 


This  formula  defines  n  implicitly  as  a  function  of  € 
and  has  been  evaluated  numerically  for  Y  *■  1.26  to  obtain 
the  tabular  relationship  used  in  DIPARDIP. 
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APPENDIX  IV 


Formulas  for  Closed  Brakes 


As  was  shown  in  the  report  "Muzzle  Brakes"  by  Hammer 
(Document  No.  7  in  the  BIBLIOGRAPHY),  the  derivation  of 
the  equations  for  closed  brake  designs  is  based  upon  find¬ 
ing  the  change  in  the  momentum  flux  thru  a  control  volume 
enclosing  the  brake. 

Regard  figure  4.1  . 


Figure  4.1  .  Schematic  Plan  View  of  Triple-Baffle,  Closed 

Ri'nlfp 


Regarding  the  passage  of  gas  toward  the  first  baffle, 
it  can  be  seen  that,  with  the  brake  flowing  full  (i.e., 
with  no  separation),  the  gas  captured  by  the  baffle  will 
have  experienced  an  approximate  area  expansion  of 

4 • 1  en  “  Sqi^+  Sn 

when  it  starts  to  be  deflected  by  the  first  baffle.  This 
gas  will  have  a  speed 

4.2  Vll|“  a°  n^€ll^  “  aonll  * 

n  is  the  speed  up  factor  developed  in  appendix  III. 

It  is  also  apparent  that  this  gas  experiences  further 
expansion  by  a  factor  S21  .  Thus,  the  resultant  total 

377 

expansion  of  the  gas  caught  and  turned  by  the  first  baffle  is 
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4.3 


€21  “  *11  S21 
511 

The  speed  of  the  gas  as  it  leaves  the  first  baffle  is 
4.4a  |v2ll“  aon21  -•  where 

b  n2i  -  n{£  21'  • 

If  the  mass  flux  thru  Sqi  is  the  same  as  that  thru 
S;q  ,  a  fraction  of  the  mass  flow  bo  equal  to 

4.5  rl  “  sll  will  be  caught  by  the  first 

S01  '+  SU 

baffle.  For  mass  conservation,  a  fraction 

4.6  K,  -  S01  will  be  passed  without  deflec- 

So;.  '+  Six 

tion  by  this  baffle. 

We  assume  that  v-q  is  nearly  paraxial  and  that  v2x 
makes  an  angle  of  <X  ^  with  respect  to  the  axis. 

For  momentum  conservation,  the  force  felt  by  the 
first  baffle  equals  the  difference  in  the  momentum  fluxes 
entering  and  leaving  it  in  the  axial  direction.  Thus, 

“0*1  |  mj  ~  | |  ccs<*t 

“oP1ac«1i  -  “oriaon21COSOCl 


fcea0rl(nil'n21COSOti)  * 

The  gas  which  passes  the  first  baffle  undeflected 
has  expanded  in  cross-sectional  area  by  a  factor  £  ]_!• 

For  full  flow  into  the  second  baffle,  the  total  area 
expansion  from  the  muzzle  is,  therefore, 

4.8  €  12  "  £  11  f  gP2_l  SL2  \ 

\  S01  /  * 

Further  expansion  of  the  deflected  gas  results  in 
an  expansion  ratio  of 

4.9  €  “  ®22  ^  12  as  the  *as  the  bafrfl 

^12 


4.7a 

b 


'  z  1 


Fzl 


or 


Fzl  " 


25 


The  velocities  of  entrance  to  and  exit  from  the  second 
baffle  corresponding  to  these  expansion  ratios  are, 
respectively, 

4.10a  V12  “  aen(€12)  "  aon12 

and 

b  v22  "  aon(€  22^  "  aon22  * 


Since  only  a  fraction  K,  of  the  total  mass  flux  enters 
the  second  baffle,  the  fraction  turned  by  this  baffle  is 


4 -11 


r«  - 


s12 


S02  +  S12 


K, 


and 


the  fraction  of  the  total  passing  undeflected  is 


4.12 


K„  -  S02 

*  5 - 


02  +  S12 


Again,  using  the  momentum  rate  of  change  principle,  we 
find  that  the  force  on  the  second  baffle  in  the  axial 
direction  is 

4,13  Fz2  “  “o  ao  r2(n12  ”  n22  COB  )  . 

By  the  same  reasoning,  the  force  on  the  third  baffle  is 


4.14 


Fz3  -  M0  a0  r3Cn13  -  n23  cos  ds  3) 


The  total  brake  force  in  the  axial  direction  is  simply 
the  sum  of  the  forces  on  each  baffle,  l.e.. 


4.15 

F  -  F  +  F  +  F 

z  zl  z2  z3 

or 

4.16 

Fz  -  «...  J  t}  <n1;)  - 

n2j  co®<^  j>  > 

where 

4.17a 

Djj  -  »<«!,) 

eii  m(sQj  4  SU  \  €ij-i 
\S0J-1 / 
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For  brakes  which  exhaust  gas  asymmetrically  with 
respect  to  the  vertical  at  zero  quadrant  elevation, 
there  is  an  additional  vertical  force  on  the  brake. 

Regard  figure  4.2  . 


Figure  4.2  .  Schematic  Front  View  of  Asymmetric,  Closed 
Brake 


with  respect  to  the  vertical.  Since  this  is  a  closed 
brake,  the  previous  analysis  which  produced  formulas  4.16 
thru  4.21  is  valid  here. 


The  downward  vertical  force  on  the  brake  is  simply 
equal  to  the  vertical  component  of  the  rate  of  change  of 
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momentum.  Since  gas  impinges  upon  each  baffle  without 
having  a  vertical  velocity  component,  the  vertical  force 
equals  the  vertical  momentum  flux.  That  is, 

4.22.  Fy  ■  5  T)  4.l»2j|co,0J  °r 

b  Fy  -  2  <j»2j  Cos/Sj  . 


It  is  convenient  to  define  several  non-dimensional 
parameters  which  are  useful  in  computing  forces  on  closed 
muzzle  brakes.  These  are: 

4.23a  X  -  F„ 

*  9 

p~ 

h  b  ■  P 

D  z  ,  the  momentum  index 

B 

c  ?\r  -  total  momentum  flux  in  axial  direction 
thru  a  control  surface  completely  en- 
closing  the  brake  _ 


A  and  b  are  alternative  ways  of  defining  a  non- 
dimensional  axial  brake  force.  ■>  r  is  a  non-dimensional 
parameter  of  interest  in  evaluating  the  strength  of  the 
induced  atmospheric  shock. 


For  closed  brakes,  Fz  is  taken  to  be  a  constant 


throughout  the  discharge  period;  hence 


4.24a  b  -  /Q  Fjj  dt  or 

"ETT" 

b  *bt-  /o  Fz  dt  * 

The  resultant  axial  impulse  felt  by  the  weapon  is 
4.25  Iz  -  B T  -  /*  Fz  dt 
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By  4.24b,  the  last  expression  may  be  written  as 

4.26  I  -  Br  (1-  b)  . 

z  1 

One  also  can  define  a  non-dimensional  vertical  brake 
force  (jj  . 

4.27  CO  -  Fy 

P 


By  integration  of  this  expression,  one  obtains  the 
resultant  vertical  impulse  Iy» 

4.28  I  -wP, 


To  find  computational  expressions  for  ,  \?  »  7\r* 

and  Ui  s  one  need  merely  make  appropriate  substitutions 
in  4.23  and  4.27  for  the  dimensional  parameters  involved. 
For  example,  from  3.21,  4.16,  and  4.23a,  one  finds 


4.29 


>-! 


r3 


Due  to  the  simplifying  assumptions  made  in  the 
derivation  of  4.29,  this  expression  is  adjusted  to  agree 
with  experiments  performed  on  closed  baffles,  by  menns 
of  a  correction  factor  1 

77 


Thus 

4.30 


N 


From  4.23b  and  the  last  expression,  the  momentum 
index  can  be  written  as 

4.31  b  “  P  1  2  ri  (n,  H-n„  .cosoC<)  . 

B  r  j  j  ^  zj  j 


The  axial  component  of  the  corrected  momentum  trans¬ 
port  thru  a  control  surface  enclosing  the  brake  can  be 
written  as 

%  *  - 
m  ctc 

T* - 

“i 


N 

i 


1  j“2Jwwo^j"“2N' 


N 

'1-2  r.) 
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With  3.21  and  4.23c  this  expression  produces 


4.32 


X 


*1 


I  rjn2jcos^j-n2N(l-|  rj) 


Finally,  using  an  empirical  correction  factor  1 


equations  3.21, 
4.33  UJ  « 


4.22b,  and  4.27  yield 
1_  |  rjn2jcos/6  j  • 
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APPENDIX  V 


Formulas  for  Open  Brakes 


For  single-baf fie,  symmetric,  open  brakes,  the 
region  between  the  muzzle  and  the  baffle  remains  rela¬ 
tively  open  for  free  expansion  of  the  gas  prior  to 
deflection.  The  baffle  is  mechanically  connected  to 
a  tube  attached  to  the  muzzle  by  supports  which  par¬ 
tially  confine  the  flow.  However,  even  without  ducting 
of  the  flow,  confinement  is  achieved  by  the  shock 
envelope  which  develops  in  flow  from  sonic  and  super¬ 
sonic  nozzles. 

llegard  figure  5.1  . 

Figure  5.1a  Schematic  Plan  View  of  a  Single-Baffle, 
Open  Brake 
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Figure  5.1b  Scl.ematic  Front  View  of  ?  Single-Baffle, 


In  the  absence  of  a  muzzle  brake,  the  free  discharge 
from  the  tube  would  cause  a  shock  envelope  to  develop  in 
the  flow.  This  has  been  sketched  in  in  figure  5.1a  .  The 
presence  of  the  brake  would  alter  the  flow  pattern  con¬ 
siderably  and  consequently  alter  any  shock  structure  within 
the  flow.  However,  it  is  not  unreasonable  to  expect  that 
the  fraction  of  the  total  mass  flux  captured  and  turned  by 
the  baffle  could  be  computed  by  assuming  that  that  portion 
is  captured  which  would  be  masked  by  a  virtual  baffle  in  a 
corresponding  position  of  the  free  stream.  With  this  as  a 
work!  fig  hypothesis,  it  is  necessary  to  be  able  to  predict 
the  free  stream  shock  envelope.  A  method  for  doing  this 
is  described  in  document  number  1  of  the  BIBLIOGRAPHY.  This 
method  was  used  for  a  sonic  nozzle  with  Hf  ■  1.26  for  several 
values  of  p0  ,  and  the  resulting  envelopes  "X  (<r  )  were 

Poo 

optimumly  fit  with  quadratic  and  cubic  equations  in  O'  . 

For  the  case  pictured  in  figure  5.1a,  that  is,  when 
ym<Q,  the  mass  fraction  captured  and  deflected  by  the 

baffle  is  given  by 

p  p 

5.1a  ri  D  Mi  *  w  yJn  -  *1  or 
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b  rj  -  1  -  iL  .  (ymSQ)  • 

yT 

However,  for  ym  >  Q  ,  the  gas  deflected  by  the  baffle 
equals  the  total  amount  minus  that  lost  thru  the  center 
and  that  which  would  be  lost  by  spilling  over  the  outer 
edge.  This  mass  is 

5.2  Mj  =*  M^_ 

?ry2 

The  fraction  captured  is 

5-3  rr  1  -  ^  t  (  J-  Q2  -  I2)/y2  ,  (ym  >  Q)  . 

n  7i 

The  momentum  flux  in  the  z-  direction  into  a  control 
surface  surrounding  the  baffle  is 

5.4  ”  v^Mi  cos  A  =  v^Mor2COS  A 

and  the  momentum  flux  in  the  z-  direction  out  of  the 
control  surface  is 

5.5  t>2  ”  Vj.foj.COS  <X  1  =  Vj^rj  cos  Ot  i 

For  momentum  conservation, 

♦  ♦ 

5.6  P1  ~  p2  a  pz  >  whi°h  from  equations  5.4  and  5.5 
may  be  written  as 

5.7  F„  -  M  r. v. (cos  A  -  cos Ot  ) 

Z  o  X  JL  1  * 

Using  the  speed-up  factor,  n,  developed  in  appendix 
III,  we  can  write 

v,  “  n  a.  ,  where 
n21  “  n<€  2V 

€  2i  "  1  +  \  and 

“»Vn  / 

A  =  0,  for  the  weapons  under  consideration. 

5.7  and  5.8a  and  d  give 

Fz  “  M0a0n21r1 (1-cos CX  1)  . 


5.8a 

b 

c 


d 

Equations 


^y2  -  ffi2  -  jt 


(yj 

wm 


Q2) 
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It  is  recognized  that  non-uniformity  in  the  density 
distribution  and  non-paraxial  flow  into  the  baffle  will 
give  rise  to  a  value  of  Fz  slightly  different  from  that 

predicted  by  equation  5.9  .  Therefore,  Fz  has  been  cor¬ 
rected  with  an  empirical  correction  factor  l/f3  .  f3  has 
been  found  to  be  approximately  1.33  .  Thus, 


5.10 

5.11 


Fz  “  a0  n21r1(l-cosa1)  or 
- - 

-  Fz  -  n2]r1(l-cosoii) 


By  4.23b  and  5.11  , 


5.12 


j?  »  P  n21r1(l-cosoti)  . 

B  13 


Since  the  open  brake  was  assumed  to  be  symmetric  with 
respect  to  a  horizontal  axis,  OJ  -  0  . 

It  should  be  noted  that  in  the  expression  for  b  , 
n21  and  rj,  both  depend  upon  the  dimensions  of  the  shock 
envelope.  Since  the  latter  changes  markedly  throughout 
the  discharge  period,  one  should  expect  to  find  some  change 
in  as  a  function  of  time  (See  graph  5.1.)  . 

Artillery  designers  are  fundamentally  concerned  with 
an  effective  momentum  index  in  order  to  compute  the  total 
impulse  which  must  be  taken  up  by  the  recoil  mechanism. 

The  engineer,  therefore,  is  not  primarily  concerned  with 
the  instantaneous  value  of  |?  .  To  find  the  effective 
value  of  which  the  engineer  uses,  we  note  that  the 
resultant  axial  impulse  is  given  by 

5.13  Iz-Bt~/*Fzdt  . 

However,  from  4.23b,  the  definition  of  b  , 

5.14  Fz  -  b  B;  hence 

ao 

5.15a  I„  “  B_.-  /  b  B  dt  or 

*  ~  0 

,Bt, 

b  I*  -  BT  -  /  bd  bt  . 
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The  effective  value  of  b  ,  ‘p  f  ,  which  the  engineer 
uses  is  obtained  from 

5.16  Iz  -  BT  (1-beff)  . 

Comparing  equations  5.15b  and  5.16,  it  is  apparent 
that  ^  should  be 

,Bt 

5.17  beff  “  . 
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Therefore*  the  downward  force  due  to  efflux  thru 
this  portion  of  the  periphery  and  its  mirror  image  with 
respect  to  the  vertical  can  be  written  as 

^2  a  n 

6.3a  Fvl  ■  2v91  K  /  .  sin  vuV-  or 


Fyl  -  2v21  k  /_jl  sin  -v-uV-  or 
Fy^  -  2v2i  K(cos"l?j-cos  where 


K  - 


For  peripheral  discharge  thru  -it  <  ’S'  <  --T&  and 

-  1 

its  mirror  image,  all  of  the  mass  is  assumed  to  be 
deflected  thru  1)3  .  Thus,  the  contribution  to  the 
downward  force  from  this  momentum  discharge  is 

6. 4  Fy2  -  2Kv21  (sinT4)(jr/2  -^1) 

Nc  flux  passes  thru  the  periphery  7^ L  -ft'  <  £_ 

and  Its  mirror  image  and,  hence,  there  is  no  vertical 
force  component  for  this  portion  of  the  periphery.  The 
total  downward  force  is,  therefore, 

6.5a  Fy  ■  2v2jK  (cosi^i-cosl^  +  Gy-  -  -zj^sint^j  0 

b  Fy  *>  M0r1v21  (cost^j-cosx^j)4^  -  T>^)sin  ^  . 

7r  L 

From  equations  4.4  and  4.27,  and  6.5b,  cO  can  be 
written  as 

6.6  co  rj_n21  J"(cos  t^-cos^2^  +  ^ - "^j)s.ln  ^ 

i  L 

For  a  multiple  baffle  brake,  this  expression  becomes 

N  i  r 

6.7  -  1_  ^  n2jrj  (cos  z  i.j“cos 

+  <*_  -'4j>sint4jj  • 

As  indicated  in  appendix  IV, 

6.8a  F„  -00$  and 


Iy  -  . 
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